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NOMENCLATURE 

empirical  constant 

empirical  constant,  sec  1 

diameter  of  cylinder 

mean  diameter  of  bed  material 

sediment  Froude  number,  /  \/ (s-l)gd 

acceleration  of  gravity 

constant  of  integration,  bed-load  movement 
constant  of  integration,  suspended-load  movement 
length  of  cylinder 

an  exponent  in  the  bed-load  transport  function 

ratio  of  ripple  amplitude  to  cylinder  diameter  at  which 

ripples  affect  the  scour 

an  exponent  in  the  suspended- load  transport  function 

O 

sediment-transport  rate,  length  3/time 
sediment-transport  rate  out  of  the  scour  hole 
ratio  of  the  specific  weight  of  the  bed  material  to 
the  specific  weight  of  the  fluid  medium 
time 

maximum  undisturbed  velocity  of  the  fluid  near  the 
bottom  above  the  boundary  layer 

vertical  distance  between  the  bottom  of  the  cylinder  and 
the  surrounding  bed  level 
initial  burial 
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NOMENCLATURE  (Continued) 

angle  of  orientation  measured  between  the  cylinder  axis 
and  the  wave  crest 
initial  orientation 

total  amplitude  of  ripple  from  crest  to  trough 
ripple  wave  length 

geometric  standard  deviation  of  particle  diameter 
angle  of  repose  of  bed  '.aterial 
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ABSTRACT 

The  settlement  of  cylindrical  mines  into  the  sea  bed  as  the  result 
of  gravity  waves  was  experimentally  studied  in  the  Georgia  Tech  Hydraulics 
Laboratory  under  general  technical  supervision  of  the  U.S.N.  Mine  Defense 
Laboratory  at  Panama  City,  Florida.  Three  s implications  were  employed  in 
the  model  tests  as  follows:  (l)  the  undisturbed  velocity  approaching  the 
mine  was  made  to  correspond  to  the  bottom  velocity  under  a  first-order 
Stokian  gravity  wave;  (2)  the  length-to-diameter  ratio  of  the  cylinder 
was  four;  and  (3)  the  bed  materials  were  limited  to  being  well-rounded  and 
uniformly  sized  having  mean  diameters  of  0.297  mm  and  O.585  mm. 

The  test  program  was  conducted  in  the  bottom  horizontal  section  of 
a  large  U  tube.  The  test  section  was  1  ft  (vertical)  by  4  ft  (horizontal) 
in  cross  section  in  which  the  water  oscillated  over  a  4-in  thick  bed  of 
sand.  Cylinder  settlement  as  a  function  of  time  was  optically  measured 
through  the  transparent  walls  of  the  test  section.  The  effects  of  initial 
submergence  and  initial  burial  were  investigated.  Settlement  was  observed 
under  waves  with  total  amplitude  of  water  motion  at  the  bed  of  2.8  ft  and 
1.0  ft  and  with  a  period  of  3-6  seconds.  Intermediate  waves  could  not  be 
utilized  in  the  model-test  program  because  of  the  formation  of  ripples 
around  the  model. 

Settlement  and  turning  was  found  to  occur  spasmodically  when  the 
support  under  the  cylinder  collapsed  as  the  surrounding  scour  hole 
deepened  and  laterally  expanded.  In  the  absence  of  ripples,  sediment 
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transport  out  of  the  scour  hole  was  observed  to  be  bed-load  movement 
with  the  lower  amplitude  waves  and  to  be  suspended- load  movement  with  the 
higher  amplitude  waves.  In  the  presence  of  ripples,  settlement  did  not 
occur  as  the  cylinder  simply  joined  the  ripple  system. 

The  data  were  analyzed  by  evaluating  the  rate  of  sediment  removal 
from  the  scour  hole.  By  consideration  of  the  forces  on  the  particles, 
by  analyzing  scour  data  of  others,  and  by  incorporating  the  model  test 
results,  rate  of  sediment  removal  functions  were  formulated  for  bed-load 
transport  and  for  suspended-load  transport.  These  transport  functions 
were  then  incorporated  into  the  equation  of  continuity  resulting  in  a 
differential  equation  of  settlement  which  was  then  integrated.  Initial 
submergence  and  initial  orientation  were  analyzed  as  being  initial  con¬ 
ditions.  In  this  manner  settlement  functions  were  formulated. 

Finally  the  settlement  functions  were  applied  to  compare  with  the 
observed  settlement  of  mines.  Data  on  mine  settlement,  bottom  velocity, 
and  bed  material  were  obtained  and  furnished  by  the  personnel  of  the 
U.S.N.  Mine  Defense  Laboratory.  Three  cases  of  mine  settlement  were 
studied  in  detail.  The  agreement  between  the  predicted  settlement  and 
the  observed  settlement  was  good  in  all  three  cases. 
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INTRODUCTION 

Mines  placed  upon  the  sea  bottom  are  known  to  bury  in  the  bottom  sands 
by  settling  into  the  encompassing  scour  hole.  The  scour  hole  results  from 
the  increased  bed-material  transport  capacity  of  the  flow  adjacent  to  the 
mine.  Mine  burial  can  be  expected  to  be  important  in  performance,  in  re¬ 
covery,  and  in  sweeping.  As  a  result  of  the  interest  and  encouragement 
from  the  Mine  Defense  Laboratory,  a  model  study  of  mine  burial  was  con¬ 
ducted  in  the  Hydraulics  Laboratory  of  the  Georgia  Institute  of  Technology 
on  Contract  NOBs -84327  for  the  U.  S.  Navy  Bureau  of  Ships,  Code  631. 

Model  studies  involving  bed-material  transport  are  characterized  by 
a  large  number  of  independent  variables  involving  flow,  geometry,  fluid, 
and  bed  material.  The  main  problem  in  an  experimental  study  is  to  elim¬ 
inate  the  independent  variables  of  negligible  influence  and  yet  to  retain 
all  those  of  major  influence.  The  writers,  by  considering  the  lift,  drag, 
fluid  inertia,  and  gravity  forces  on  particles  resting  on  the  surface  of 
the  bed,  concluded  that  a  sediment  Froude  number  is  the  primary  similarity 
parameter  involving  sediment  properties  (Appendix  A).  The  modelling 
relation  for  the  cylindrical  mine  on  a  movable  bed  is 
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cylinder  diameter; 
cylinder  length; 

initial  orientation  of  cylinder; 

ratio  of  the  bed -material  density  to  the  fluid  density; 
mean  diameter  of  the  bed  particles;  and 
geometric  standard  deviation  of  particle  diameter 


The  model  tests  were  performed  using  three  simplifications.  First 
the  undisturbed  velocity  approaching  the  mine  was  made  to  be  simple 
harmonic  corresponding  to  the  bottom  velocity  under  a  first-order  gravity 
wave.  Second,  the  model  was  a  right  circular  cylinder  with  a  length-to- 
diameter  ratio  of  four.  Third  only  uniformly-sized  bed  materials  were 
used  thereby  eliminating  the  variable  c  from  the  experimental  program. 

The  model  tests  were  performed  in  the  horizontal  bottom  leg  of  a 
large  U-tube.  The  amplitude  of  oscillation  within  the  U-tube  could  be 
controlled  but  the  frequency  was  constant.  The  cylindrical  model  was 
placed  on  the  bed  at  various  angles  of  orientation,  ,  and  with  various 

initial  submergence  levels,  ygl/D.  A  record  of  burial  was  obtained  by 
means  of  a  cathetometer  during  the  majority  of  mans  and  by  means  of  l6-mm 
motion  pictures  during  two  runs .  Four  of  uhe  independent  variables  of 
equation  (l),  namely,  U^t/D,  y^/n,  and  oC  q  were  varied  as  desired 
during  the  experiments . 

The  sediment  Froude  number,  Um/  s/T^-l)  gd,  could  not  be  varied 
throughout  the  desired  range  even  though  the  velocity,  U  ,  and  the  sedi¬ 
ment  diameter,  d,  could  be  varied.  With  intermediate  values  of  the 
sediment  Froude  number,  ripples  formed  on  the  bed.  The  appearance  of 
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ripples  on  the  bed  of  the  same  order  of  magnitude  as  the  model  precluded 
performing  model  tests  with  a  rippled  bed.  Early  in  the  experimental 
program  the  discovery  was  niaue  that  ripples  would  form  at  conditions  less 
than  the  incipient -motion  condition  with  the  result  that  the  flow  dis¬ 
turbance  created  by  tne  model  would  initiate  ripples  which  would  then 
spread  over  the  bed.  The  unexpected  appearance  of  ripples  severely 
limited  the  range  through  which  the  sediment  Froude  number  could  be 
varied.  Consequently  tests  were  made  at  conditions  of  low  wave  ampli¬ 
tude  and  at  conditions  of  high  wave  amplitude  in  order  to  avoid  ripples . 

A  detailed  analysis  of  scour  was  executed  (Appendix  A)  in  order  to  inter¬ 
polate  the  test  results  for  the  intermediate  wave  amplitudes . 


CONFIDENTIAL 


3 


CONFIDENTIAL 


EXPERIMENTAL  PROGRAM 

In  order  to  study  scour  on  the  sea  bed  resulting  from  wave  action, 
the  decision  was  made  to  model  only  the  mass  of  the  water  adjacent  to  the 
bed.  The  water  motion  at  a  fixed  point  close  to  the  bed  under  a  first- 
order  Stokian  wave  is  simple  harmonic  and  is  parallel  to  the  bed.  A  large 
U-tube  with  forced  oscillation  of  the  water  was  designed  in  order  to  model 
the  water  motion  under  a  wave. 

U-Tube 

Description  -  The  description  of  this  large  U-tube  is  facilitated  by 
referring  to  Figure  1.  The  vertical  legs  of  the  U-tube  are  in  two  rec¬ 
tangular  steel  tanks  (A)  at  the  ends  of  the  horizontal  leg  which  is  the 
test  section  (b).  Forced  oscillation  of  the  water  mass  is  achieved  by 
blowing  air  into  the  West  vertical  leg  as  the  water  surface  is  falling 
and  then  exhausting  this  air  as  the  water  surface  is  rising. 

The  vertical  legs  of  the  U-tube  are  formed  within  the  rectangular 
steel  tanks  which  are  3  ft  by  4  ft  in  cross  section,  by  streamlined 
inserts  (C) .  The  water  passage  in  each  vertical  leg  is  1  ft  by  4  ft  in 
cross  section  inasmuch  as  the  water  surface  is  never  allowed  to  fall  to 
the  curved  section  of  the  upper  insert  (C).  In  all  tests  the  equilibrium 
water  level  was  established  48-1/2  in  above  the  top  of  the  test  section  (b) . 

The  horizontal  leg  of  the  U-tube  is  the  test  section  which  is  1  ft 
(vertical)  by  4  ft  (horizontal)  in  cross  section  and  which  is  10  ft  long. 

The  central  portion  of  the  floor  is  depressed  in  order  to  form  a  container 
for  the  erodible  bed  material.  The  erodible  bed  (D)  is  6  ft  long  by  4  ft 
wide  by  4  in  deep.  The  walls  of  the  test  section  are  fabricated  of  1/2-in 
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Figure  1.  Side  Election  and.  Cross  Section  of  Tube 
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clear  plastic  and  are  framed  on  the  exterior  with  steel  angles  and  channels . 
The  test  section  rests  upon  three  prefabricated  steel  trusses  which  span 
from  steel  tank  (A)  to  steel  tank.  A  3-ft  square  flush-mounted  door  is 
located  in  the  center  of  the  roof  in  order  to  he  able  to  place  the  bed 
material  and  models . 

The  water  in  the  U-tube  is  made  to  oscillate  at  the  resonant  frequency. 
The  output  of  a  centrifugal  blower  is  discharged  continuously  into  the  air 
space  above  the  water  surface  of  the  West  vertical  leg.  Two  7-in 
diameter,  pneumatically  powered,  exhaust  valves  (E)  in  the  top  of  West 
vertical  leg  are  opened  upon  receipt  of  the  signal  that  the  minimum  water 
level  in  this  leg  is  attained. 

The  feedback  mechanism  by  which  the  exhaust  valve  is  sequence-operated 
at  the  resonant  frequency  is  as  follows.  The  float  (P)  in  the  East  vertical 
leg  is  attached  by  a  light  flexible  cable  to  a  steel  rod  (G)  which  moves 
vertically  past  the  direction-sensing  switch  (h).  The  direction-sensing 
switch  (H)  is  a  lever-operated  microswitch.  A  permanent  magnet  on  the  end 
of  the  microswitch  operating  lever  is  in  contact  with  the  steel  rod  (g) 
which,  in  turn,  follows  the  motion  of  the  float  (F).  Whenever  the  steel 
rod  (g)  is  rising  the  switch  (h)  is  closed  and  whenever  the  steel  rod  (g) 
is  falling  the  switch  (H)  is  open.  When  the  steel  rod  (g)  changes 
direction  and  starts  to  rise,  a  circuit  is  closed  which,  in  turn,  actuates 
a  single-cycle  timer.  This  timer  makes  one  revolution  in  2  seconds  and 
then  stops.  A  second  microswitch  is  contained  within  the  timer.  3y  means 
of  an  adjustable  cam  this  second  microswitch  can  be  made  to  open  or  close 
at  any  time  within  the  two-second  interval.  Solenoid  valves  which  operate 
the  pneumatic  pistons  on  the  exhaust  valves  are  in  the  circuit  with  the 
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timer  microswitch.  The  timer  microswitch  is  set  such  that  the  exhaust 
valves  open  when  the  timer  starts  and  such  that  the  exhaust  valves 
remain  open  for  a  half  period.  The  feedback  mechanism  described  above 
insures  that  the  water  is  oscillated  at  resonant  frequency  with  the 
result  that  the  frequency  of  oscillation  cannot  be  controlled. 

An  electrically  operated  counter  is  placed  in  the  circuit  containing 
the  direction-sensing  switch  (h)  for  the  purpose  of  indicating  the  number 
of  oscillations  or  waves  from  the  beginning  of  each  run. 

The  amplitude  of  the  oscillation  can  be  controlled  by  positioning 
the  cone  valve  placed  on  the  inlet  to  the  centrifugal  blower.  For  the 
determination  of  amplitude  a  scale  is  fixed  parallel  to  the  steel  rod  (g). 
A  pointer  on  the  steel  rod  passes  over  the  face  of  the  fixed  scale. 

A  blow-down  system  and  a  quick-opening  flap  valve  (l)  were  installed 
in  the  East  vertical  leg  in  order  to  eliminate  the  transients  associated 
with  starting  the  desired  oscillation.  A  solenoid  valve  in  the  high-pres¬ 
sure  air  line  leading  to  the  top  of  the  East  vertical  leg  permitted  the 
water  surface  to  be  depressed  to  the  lowest  elevation  of  the  final  oscil¬ 
latory  motion.  The  flap  valve  is  fitted  with  a  quick-release  mechanism 
in  order  to  be  able  to  begin  the  oscillation  when  the  water  surface  was 
blown  down  to  the  desired  value.  Upon  release  of  the  flap  valve ,  the 
feedback  mechanism  associated  with  the  direction-sensing  switch  is 
activated  with  the  steady  cyclic  motion  being  maintained  thereafter. 

The  photograph,  Figure  2,  shows  the  South  wall  of  the  test  section. 

Calibration  -  The  operating  characteristics  of  the  U-tube  and  the 
flow  characteristics  within  the  test  section  have  been  reported  in  a 
separate  report  (l).  The  principal  conclusions  were  (l)  that  the  motion 
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Figure  2 .  Photograph  of  U-Tu'oe . 
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is  simple  harmonic  within  the  test  section,  (2)  that  the  fluid -particle 
amplitude  within  the  test  section  is  equal  to  the  float  amplitude ,  and 
(3)  that  the  velocity  distribution  is  essentially  uniform  within  the 
test  section. 

Models 

The  right-circular  cylinder  models  are  fabricated  from  aluminum  with 
diameters  of  0.501  in,  1.002  in,  1-702  in,  and  3 • ^50  in  and  with  a  length - 
to -diameter  ratio  of  four.  Perpendicular  diameter  markings  on  the  ends 
connected  by  four  axial  markings  on  the  curved  surface,  Figure  3>  arc 
aids  in  visual  determination  of  displacement. 

Bed  Material 

The  bed  materials  are  either  glass  beads,  catalog  number  090,  obtained 
from  the  Minnesota  Mining  and  Manufacturing  Company  or  "Flint  Shot"  Ottawa 
sand  obtained  from  the  Ottawa  Silica  Company .  The  principal  advantage  of 
these  materials  is  their  uniformity  in  regard  to  diameter.  Size  character¬ 
istics  were  determined  by  sieving.  Specific  gravity  was  determined  by  a 
standard  pycnometer  technique.  Angle  of  repose  was  determined  by  repeated 


measurements 

of  the  slope  of 

a  submerged  pile  of  the 

sand  which  was  formed 

by  pouring  between  parallel 

glass  plates. 

The  characteristics  are  tabu- 

lated  below 

Material 

Mean  Diameter 
d  (mm) 

Geometric 
Std.  Dev. 

°”p;d 

Specific 

Gravity 

s 

Angle  of  Repose 
(deg) 

Glass  Beads 

0.297 

1.06 

2.47 

2h 

Ottawa  Sand 

0.585 

l.l6 

2.62 

32.5 
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Experimental  Procedure 

Prior  to  making  a  run  the  following  operations  were  accomplished. 
First,  the  sand  bed  was  leveled  with  a  2-in  by  2-in  wooden  screed  which 
spanned  the  six-foot  length  of  the  bed.  The  sand  was  kept  saturated 
during  the  leveling.  An  external  vibrator  was  moved  over  the  exterior 
of  the  test  section  to  facilitate  compaction  of  the  bed  and  to  facilitate 
removal  of  trapped  air.  When  the  bed  was  level,  water  was  admitted  until 
the  water  surface  was  approximately  6  in  above  the  bed  surface.  A  string 
was  then  stretched  across  the  inside  of  the  tunnel  just  above  the  water 
surface  as  an  aid  in  aligning  the  model.  A  model  was  then  gently  lowered 
to  the  bed  matching  the  string  and  the  axial  marks  on  the  model  in  order 
to  attain  the  desired  orientation.  Then  the  door  in  the  roof  was  re¬ 
placed  and  rebolted  in  position.  Next  the  U-tube  was  filled  with  water 
to  equilibrium  level. 

After  the  model  had  been  positioned  in  the  test  section,  a  cathe- 
tometer  with  a  20X  telescope  was  placed  adjacent  to  the  South  wall  of 
the  test  section.  Elevations  of  the  top  of  the  sand  bed  and  top  of  the 
model  were  then  determined  and  recorded. 

Four  observers  were  required  during  a  run.  The  functions  of  the 
first  observer  were  (l)  to  blow  down  the  East  leg  to  a  predetermined 
level,  (2)  to  release  the  flap  valve  at  the  top  of  the  East  leg,  (3) 
to  flash  a  300-watt  lamp  at  selected  intervals  (generally  5  cycle 
intervals),  (4)  to  read  and  record  amplitude,  and  (5)  to  read  and  re¬ 
cord  elapsed  time  for  a  given  number  of  cycles.  The  second  observer 
manned  the  cathetometer  telescope  maintaining  the  cross  hairs  on  the 
top  of  the  model.  The  third  observer  read  and  recorded  the  cathetometer 
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reading  when  the  first  observer  flashed  the  lamp.  The  fourth  observer 
was  stationed  on  top  of  the  test  section  with  a  large  protractor .  The 
fourth  observer  read  and  recorded  the  model  orientation  as  the  signal 
lamp  was  flashed. 

Water  temperature  was  read  and  recorded  either  just  before  or  just 
after  a  run. 

The  procedure  was  altered  during  Model  tests  8  and  43  when  the 
model  settlement  was  recorded  by  a  l6-mm  motion  picture  camera.  During 
these  two  runs  a  running  clock  was  placed  in  the  field  of  view  of  the 
camera.  Burial-time  data  were  determined  by  scaling  from  projected 
images  of  the  l6-mm  colored-film  strip. 

The  fourth  observer  was  not  required  during  the  runs  where  the  model 
axis  was  initially  parallel  to  the  wave  crest,  that  is,  for  runs  in  which 

=  o. 

o 

Results 

All  model-test,  data  are  tabulated  and  enclosed  in  the  colored  sheets, 
pages  13-48,  inclusive. 
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MODEL  TEST  5 

BED  MATERIAL — GLASS  BEADS  TEMPERATURE=65.0°  F 


TOTAL  AMPLITUDE*32.60  IN  D=1.702  IN  =15° 

© 


CYCLES 

Y/D 

SOUTH 

CYCLES 

Y/D 

SOUTH 

0 

0.127 

19 

0.717 

4 

0.430 

24 

0.836 

9 

0.594 

29 

0.926 

14 

0.717 

34 

0.970 
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MODEL  TEST  6 


BED  MATERIAL*— GLASS  BEADS  TEMPERATURE=64»Q°  F 

TOTAL  AMPLITUDE*33fi30  IN  D=1.702  IN  <*o  *30° 


CYCLES 

Y/D 

SOUTH 

CYCLES 

Y/D 

SOUTH 

0 

0*632 

14 

0*836 

4 

6*490 

19 

0*  $89 

9 

0*662 

24 

0.975 
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MODEL  TEST  7 


BED  MATERIAL*— GLASS  BEADS  TEMPERATURE=64.(f  F 

TOTAL  AMPLITUDES  3  •  55  IN  D=1.702  IN  c*=6 <f 


CYCLES 

Y/D 

SOUTH 

CYCLES 

Y/D 

SOUTH 

0 

O«058 

14 

0.773 

4 

0*548 

19 

0*850 

9 

0.620 

24 

0*940 

29 

0*954 
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BED 
TOTAL  t 
CYCLES 


0*000 
0#  }22 
0.268 
01346 
0.56$ 

o.ese 

0*805 

1*032 

i*2it 

1.317 

1*545 

1*816 

2*040 

2.440 

2.560 

3.067 

3.580 


MODEL  TEST  8 


MATERIAL— GLASS  BEADS 


TEMPERATURE=67. 0°  F 


IMPLITUDE«34»20  IN 


D=1.702  IN 


*0° 


Y/D 

CYCLES 

Y/D 

SOUTH 

SOUTH 

0.040 

4.070 

0  j  490 

01049 

4.575 

0.507 

0*064 

5.09§ 

0.507 

0.067 

5.660 

0*307 

0.131 

6.005 

0*615 

0.119 

7.085 

0.558 

G.1J7 

8.060 

0.643 

0.303 

8.700 

0.677 

0*396 

10. 070 

0.694 

0.345 

12.050 

0.728 

0*363 

14.000 

0.762 

0*363 

16.010 

0.813 

0.439 

18.530 

0.813 

0.465 

20.570 

0.848 

0.473 

23.480 

0.915 

0.490 

26.732 

0.949 

0.490 

27.400 

0.966 

49.000 

1.197 
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MODEL  TEST  9 


MATERIAL— GLASS  BEADS  TEMPERATURE-67. £  F 


AMPLITUDE =34. 10  IN 


D=3*45  IN 


«  ■  14 


Y/D 

SOllTH 


CYCLES  Y/D 

SOUTH 


0*062 

0.142 

0.171 

0*251 

0*411 

0*463 

0.621 

0*643 


39 

0*713 

44 

0.760 

49 

0*800 

54 

0*819 

59 

0*848 

64 

0*9?2 

69 

0*904 

74 

0*910 

sp 

0*$53 
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MODEL  TEST  10 


BED  MATERIAL — GLASS  BEADS 
TOTAL  AMPLITUDE-34.00  IN  D-3.45 


YCkES 

Y/D 

SOUTH 

0 

0*078 

4 

0*257 

9 

9*428 

14 

0*519 

19 

0*632 

24 

0*664 

29 

0*691 

34 

0*728 

TEMPERATURE-69. 0°  F 
IN  -60° 

O 


CYCLES 

Y/D 

SOUTH 

39 

0*761 

44 

0*860 

49 

0*857 

54 

0*964 

59 

0*947 

64 

0*973 

69 

1*661 

74 

1*021 

91 

i*i£o 
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model  test  11 


BfD  MATERIAL — GLASS  BEADS 

TOTAL  AMPLITUDES. 70  IN  | 

CYCLES  Y/D 

SOUTH 


0 

9 

14 

19 

24 

29 

34 

39 


0*077 

0*312 

0*420 

0*477 

0*566 

0*587 

0*679 

0.700 


TEMPERATURE-70*  0°  F 


*<o“30' 


CYCLES 


44 

49 

54 

59 

64 

69 

74 

n 


Y/D 

SOUTH 

0*731 

0.760 

0.793 

0*643 

0.882 

0.913 

0.929 

0.940 
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MODEL  TEST  12 


BED 

MATERIAL— GLASS  BEADS 

TEMPERATURE* 

67.0°  F 

TOTAL 

AMPLITUDES. 80  IN  0=3.45 

IN  0^=15° 

CYCLES 

Y/D 

CYCLES 

Y/D 

SOMTM 

SOUTH 

0 

0.082 

44 

0.816 

4 

0.312 

49 

0.824 

9 

0*387 

54 

0.830 

14 

0.477 

59 

0.855 

19 

0.571 

64 

0.&86 

24 

0.607 

69 

0.924 

29 

0.663 

74 

0.961 

34 

0.730 

79 

0.981 

3$ 

0.798 

88 

1.000 
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MODEL  TEST  13 


BED  MATER JAL--GLASS  BEADS  TEMPERATURE*70.0°  F 

TOTAL  AMPLITUDE*33.60  IN  D*3«45  IN  o(  •q<3 


CYCLES 


0 

5 

10 

15 

20 

25 

30 

35 


Y/D 

SOUTH 

0.078 

0t524 

0*369 

0*474 

0*540 

0.528 

0*618 

0*711 


CYCLES 


40 

45 

50 

55 

60 

65 

70 

75 

90 


Y/D 

SOUTH 

0*732 

0*759 

6*799 

0*851 

0*892 

0»9?7 

0*956 

0*987 

1*007 
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MODEL  TEST  14 


BED 

MATERIAL— GLASS  BEADS 

TEMPERATURE’ 

TOTAL 

AMPLITUDE *33 *70  IN  D*3.45 

IN  *9{f 

CYCLES 

Y/D 

CYCLES 

WEST 

0 

0.078 

70 

5 

0*172 

75 

10 

0.290 

80 

1$ 

0*245 

90 

20 

0*281 

95 

25 

0*308 

109 

30 

0*324 

110 

35 

0.335 

115 

40 

0. 33$ 

120 

50 

0.346 

130 

55 

0.351 

145 

60 

0.357 

155 

45 

0*364 

10 

64*2°  F 


Y/D 

WEST 

0.375 
0.384 
0.369 
0 . 4Q7 
0.413 
0.431 
0.439 
0.445 
0.455 
0.485 
0.491 
0.50* 
0.520 
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MODEL  TEST  15 


BED  MATERIAL— GLASS  BEADS  TEMPERATURE=65 . 2°  F 


TOTAL 

AMPLITUDE*34*00  IN 

D=3.45  IN  *75° 

CYCLES 

Y/D 

CYCLES 

Y/D 

SOUTH 

SOUTH 

0 

0*089 

30 

0*663 

5 

0*166 

35 

0*729 

10 

0.291 

40 

0*774 

15 

Q.3$7 

45 

0*891 

20 

0.493 

SO 

0.893 

25 

0.610 

55 

0«942 
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MODEL  TEST  16 


BED 

MATER  I AL-aGLASS  BEADS 

TEMPERATURE-66. 0°  F 

TOTAL 

AMPLI TUDE*33 *80  IN  D=3*45 

IN  <*„ 

*80° 

cycles 

Y/D 

CYCLES 

Y/D 

EAST 

EAST 

0 

0  *  065 

45 

0*712 

5 

0*226 

50 

0*747 

10 

0*334 

55 

0*772 

15 

0.397 

60 

0*777 

20 

0*436 

65 

0*782 

?5 

0.479 

70 

0*800 

30 

0*373 

73 

0.794 

35 

0*616 

80 

0*  800 

4(3 

0*658 

95 

0.917 
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MODEL  TEST  17 


BED  MATERIAL— GLASS  BEADS 


total 

AMPLITUD£*34.40  IN 

cycles 

Y/D 

EAST 

0 

0*145 

5 

0.633 

TEMPERATURE^?. 5  F 
0=1*702  IN  «75° 


CYCLES 

Y/D 

EAST 

10 

0*747 

15 

0*9?1 

CONFIDENTIAL 


25 


CONFIDENTIAL 


MODEL  TEST  18 


BED  MATERIAL— GLASS  BEADS 


TOTAL 

AMPLITUDE -34. 05  IN 

CYCLES 

Y/D 

WEST 

0 

0.100 

5 

0.476 

10 

0.586 

15 

0.638 

TEMPERATURE-66. 0°  F 
=  1.702  IN  =89.5° 


CYCLES 

Y/D 

WEST 

20 

0.702 

25 

0.791 

30 

0*811 

35 

0*837 
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MODEL  TEST  24 


BED  MATERIAL— -OTTAWA  SAND  TEMPERATURE-74. 3°  F 


TOTAL  AMPLITUDE-12.15  IN  D-0.501  IN 


CYCLES  Y/D 


SOUTH 

NORTH 

0 

0.075 

0,075 

5 

0*378 

10 

0.551 

15 

0,567 

20 

0.575 

25 

0*583 

30 

0.630 

35 

0.818 

0.433 

40 

0.685 

45 

0.6?3 

56 

0.709 

55 

0.709 

60 

0.717 

55 

0.748 

ip 

0.748 

CYCLES 

Y/D 

SOUTH 

NORTH 

75 

0.551 

80 

0.764 

85 

0.787 

90 

0.787 

95 

0.803 

100 

0.808 

03 

0.827 

115 

0.843 

130 

~  i 

0.630 

140 

0.827' 

145 

0.906 

150 

0.921 

160 

0.921 

170 

0*992 

180 

0.606 

190 

l.OOQ 
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MODEL  TEST  25 


BED  MATERIAL— OTTAWA  SAND 


total 

amplitude 

=12.09 

CYCLES 

y/d 

SOUTH 

NORTH 

0 

0.032 

0.032 

5 

0.350 

15 

0*216 

20 

0.555 

25 

0.087 

30 

0.507 

35* 

0.567 

40 

0.106 

49 

0.606 

50 

0.854 

55 

0.036 

60 

0.874 

65 

0.874 

70 

0.874 

75 

0.098 

TEMPERATURE *73  .  5* F 
D=0.501  IN  c* 


CYCLES 

SOUTH 

Y/D 

NORTH 

80 

90 

0.831 

0.240 

95 

160 

105 

0.854 

0.240 
0.27  2 

110 

0.299 

120 

1  DR 

ft  ft  *  A 

0.331 

130 

0*370 

135 

140 

0.890 

0.386 

145 

150 

155 

160 

170 

0.894 

0.380 

0.417 

0.417 

0.417 

0.417 
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MODEL  TEST  26 


BED  MATERIAL— OTTAWA  SAND 


total 

AMPLITUDE 

*11.90 

CYCLES 

Y/D 

SOUTH 

NORTH 

0 

0.004 

5 

0.402 

10 

* 

0.295 

1? 

0.406 

20 

‘  ‘  '  ** 

0.248 

25 

0.630 

30 

0.197 

35 

0.697 

40 

•  -  1 

0*197 

45 

0.760 

50 

0.197 

55 

0.764 

60 

0*764 

65 

0.764 

70 

1 

0.331 

7? 

0.705 

TEMPERATURE=77.0  F 
D=0»501  IN  ex'  =o* 


CYCLES 

85 

SOUTH 

0.894 

Y/D 

NORTH 

90 

0.264 

95 

0.886 

160 

no 

0.886 

0.886 

115 

120 

0.870 

0.327 

125 

0.870 

130 

135 

140 

145 

150 

0.941 

0.402 

0.413 

0.413 

0.413 

160 

165 

0.953 

0.429 

170 

0.953 

0.429 
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MODEL  TEST  27 


BED  MATERIAL*— OTTAWA  SAND  TEMPERATURE=78«4e>  F 


total 

AMPLITUDE 

=12.14  IN 

D=0»501  IN  <V 

■o’5 

CYCLES 

Y/D 

CYCLES 

Y/D 

SOUTH 

NORTH 

SOUTH  NORTH 

0 

0.077 

0.077 

180 

0.416 

5 

0.272 

0.272 

185 

0.496 

10 

0.272 

190 

0.406 

15 

0.350 

195 

0*406 

20 

•  *; 

0.280 

200 

0.406 

25 

0.496 

205 

6.429 

30 

0.515 

210 

0.429 

35 

■* 

0.260 

215 

0.437 

40 

0.606 

220 

0.437 

45 

0.622 

2*5 

6.437 

50 

0.681 

230 

0.437 

55 

0.295 

23? 

0.437 

60 

0.70? 

240 

0.437 

65 

0.343 

245 

0.437 

70 

0.68? 

250 

0.437 

75 

" 

0.3?2 

255 

6.437 

80 

0.669 

260 

6.437 

8? 

0.382 

265 

0.437 

90 

0-823 

270 

0.437 

95 

0.315 

275 

0.437 

100 

0.906 

280 

0.390 

105 

0.299 

285 

0*390 

no 

0.827 

290 

6.390 

115 

0*311 

297 

0*453 

120 

0.835 

365 

0.539 

125 

0.327 

310 

0.575 

130 

0*980 

315 

6.594 

135 

0.350 

320 

p.594 

140 

0.988 

‘  '  " 

325 

0.594 

145 

0.354 

330 

0.657 

150 

0.996 

335 

0.669 

155 

0.336 

340 

0.693 

160 

0.980 

345 

0.748 

165 

0.406 

350 

0*823 

170 

0.969 

355 

0*823 

175 

0.437 

358 

1.449  0.776 
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MODEL  TEST  28 


BED  MATERIAL — OTTAWA  SAND  TEMPERATURE  =  79. 5°  F 


TOTAL 

AMPLITUDE 

“12.23  IN 

D=0. 501  IN 

O 

*0 

CYCLES 

Y/D 

CYCLES 

Y/D 

SOUTH 

NORTH 

SOUTH 

NORTH 

0 

0*274 

0.274 

150 

0.634 

2 

0.378 

155 

1.067 

5 

0*465 

160 

0*665 

10 

0.603 

165 

1.028 

15 

0.701 

170 

0.720 

20 

0.717 

175 

1.020 

25 

0*736 

180 

0.728 

30 

0.213 

185 

1.035 

?s 

0.791 

190 

0.736 

40 

0.272 

195 

1.051 

45 

0.776 

200 

0.756 

50 

0.280 

205 

1.051 

S5 

0.791 

210 

•  1 

0.748 

65 

1.012 

220 

0.756 

70 

0.228 

225 

1.07? 

73 

0.998 

230 

0.748 

80 

0.287 

235 

1.Q79 

•  •  , 

85 

0.996 

245 

0.756 

90 

0.413 

256 

1.06? 

99 

1.087 

255 

0.736 

100 

0.413 

265 

0.736 

106 

1.094 

270  ' 

1.150 

■  •• 

110 

0.492 

279 

0.720 

115 

1.047 

280 

1.173 

120 

0.543 

285 

0.882 

125 

1.122 

290 

0.831 

130 

0*531 

295 

0*839 

135 

1.138 

300 

0.839 

140 

0.567 

330 

0.839 

145 

1.106 

355 

1.590 

0.83? 
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MODEL  TEST  29 


BED  MATERIAL — OTTAWA  SAND 
TOTAL  AMPLITUDES. 16  IN 


CYCLES  Y/D 


SOUTH 

NORTH 

0 

0.360 

0.360 

2 

0.472 

0.472 

5 

0.472 

10 

0.474 

0.474 

15 

0.492 

20 

0.460 

25 

0.528 

10 

0.543 

35 

0.559 

46 

0.555 

0.555 

45 

0.559 

0.559 

50 

0.596 

55 

0.587 

6  0 

0.665 

65 

0.587 

70 

0.701 

- 

76 

0.559 

80 

0.756 

85 

0.555 

95 

0.788 

100 

0.508 

105 

0.906 

' 

110 

0.472 

115 

0.976 

126 

0.472 

125 

0.969 

130 

0.465 

135 

0.967 

140 

0.474 

150 

1.118! 

TEMPERATURE =78 * 5  F 
D=0»  501  IN  o<  =o° 

O 

CYCLES  Y/D 


SOUTH 

NORTH 

155 

0*600 

160 

1.122 

165 

0.504 

170 

1.051 

173 

0.543 

180 

1.102 

185 

0.583 

190 

1.071 

1 

195 

0.591 

200 

1.105 

205 

0.610 

210 

0.610 

215 

0.610 

225 

0.646 

230 

0.646 

235 

0*646 

240 

0.689 

245 

0.689 

250 

0.689 

255 

0.713 

275 

0.713 

300 

0.713 

305 

0.724 

310 

0.724 

315 

0.724 

320 

0.748 

330 

0.748 

355 

0.748 

370 

0.748 

375 

1.201 

0.740 
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MODEL  TEST  30 


BSD  MATERIAL— OTTAWA  SAND  TEMPERATURE=78. 0°  F 


TOTAL 

AMPLITUDE 

=12.00  IN 

D  =  0»  501  IN  o^o 

“0° 

CYCLES 

Y/D 

CYCLES 

Y/D 

SOUTH 

NORTH 

SOUTH 

NORTH 

0 

0»551 

0.551 

180 

0.906 

2 

0.555 

0.555 

183 

0.906 

5 

0*555 

0.555 

190 

0.906 

10 

0.559 

0.559 

195 

0.906 

15 

0.559 

0.559 

200 

0.906 

20 

0.559 

0.559 

205 

0 .906 

25 

0.559 

0.559 

210 

' 

0.630 

30 

0.559 

0.559 

220 

1.067 

36 

0.559 

0.559 

225 

0.559 

40 

0.559 

0.359 

230 

1.079 

45 

0.543 

235 

0.579 

50 

0.614 

240 

0.579 

55 

0.575 

245 

Qe  5  04 

60 

0.618 

0.618 

250 

0.504 

65 

0.661 

255 

0.504 

70 

0.689 

0.689 

260 

0.520 

75 

0.717 

0.717 

265 

0.520 

80 

0.764 

270 

0.520 

85 

0.701 

275 

0.520 

90 

0.764 

280 

0.534 

95 

0.764 

285 

0.540 

100 

0.772 

290 

0.574 

105 

0.697 

295 

0.614 

110 

0.779 

300 

0*618 

115 

0.811 

305 

0.624 

120 

0.689 

310 

0.701 

125 

0.874 

315 

0.760 

130 

0.890 

320 

0.764 

135 

0.677 

325 

0.784 

140 

0.898 

330 

0.795 

145 

0.638 

335 

0.811 

150 

0.902 

340 

0.819 

155 

0.634 

345 

0.827 

160 

0.902 

350 

0.854 

165 

0.902 

355 

0.890 

170 

0.902 

360 

0.937 

175 

0.630 

365 

1.405 

0.961 
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MODEL  TEST  31 


BED  MATERIAL— -OTTAWA  SAND 


TEMPERATURE*77.5°  F 


TOTAL  AMPLITUDE =11, 9 3  IN  D*1.00  IN 


CYCLES  Y/D 

SOUTH  NORTH 

0  0,060  0.060 

2  0.123  0.123 

5  0.210  0.210 

10  0.273  0.273 

RIPPLES  STARTED  FORMING, 


CYCLES 

Y/D 

SOUTH 

NORTH 

15 

0.279 

0.279 

20 

0.264 

25 

0.278 

30 

0.284 

TEST  DISCONTINUED 
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MODEL  TEST  32 


BED 

MATERIAL — OTTAWA  SAND 

TEMPERATURE=76. 

o 

5  F 

TOTAL 

AMPLITUDE*^. 25  IN  D=0.501  IN 

L/D  *  8 

ex’ 

o 

*0° 

CYCLES 

Y/D 

CYCLES 

Y/D 

SOUTH  NORTH 

SOUTH 

NORTH 

0 

0.063 

160 

0.972 

2 

0.150 

165 

0.410 

5 

0.417 

170 

0.997 

10 

0.343 

175 

0.968 

15 

0.417 

180 

0.441 

20 

0.394 

185 

''.961 

25 

0.402 

190 

0.961 

30 

0.402 

195 

0.977 

35 

Q.500 

2Q0 

0.980 

46 

0.390 

205 

0*969 

45 

0.587 

210 

0.512 

50 

0.515 

319 

1  i  QQO 

55 

0*327 

220 

0.499 

60 

0.701 

225 

0.563 

65 

0.322 

230 

0*563 

70 

0.715 

235 

0.563 

75 

0.815 

240 

0.563 

80 

0.325 

245 

0.510 

85 

0.803 

255 

1.016 

90 

0.842 

260 

0.504 

95 

0.394 

265 

0.504 

100 

0.852 

270 

1.016 

105 

0.852 

275 

1*016 

110 

Q.  363 

280 

0.512 

115 

0.890 

285 

0.512 

126  1 

0.881 

290 

0.512 

125 

0.386 

295 

0*512 

130 

0.922 

300 

0.512 

135 

0.922 

305 

0.543 

140 

0.905 

310 

0.543 

145 

0.386 

315 

0.543 

150 

0.933 

340 

0.543 
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MODEL  TEST  33 


BED  MATERIAL — OTTAWA  SAND  TEMPERATURE=76« 5° F 

L/D  -  12 


TOTAL  AMPLITUDES. 23  IN  D=0.501  IN 


CYCLES 

Y/D 

CYCLES 

SOUTH 

NORTH 

SOUTH 

0 

0.145 

235 

0.732 

2 

0.241 

240 

0.673 

5 

0.421 

255 

0.685 

10 

0.437 

260 

0.733 

15 

0.437 

265 

20 

0.437 

270 

0.764 

25 

0.437 

275 

0.764 

30 

0.437 

280 

0.740 

$5 

0.437 

285 

40 

0.437 

290 

/ 1717 

45 

0.437 

295 

733 

50 

0.437 

300 

0.799 

55 

0.437 

305 

60 

0*437 

310 

0.748 

65 

0.437 

315 

0.748 

70 

0.437 

320 

0.748 

75 

0.437 

325 

0.748 

80 

0.437 

330 

0.819 

85 

0.495 

335 

90 

0.520 

340 

0.827 

95 

0.599 

345 

0.827 

100 

0.500 

350 

0.795 

110 

0.481 

355 

0.795 

120 

0.602 

360 

0. 850 

125 

0*402 

365 

130 

0.527 

370 

0.847 

135 

0.591 

375 

0.847 

140 

0.527 

380 

0.772 

145 

0.504 

385 

0.77? 

150 

0*650 

395 

0.906 

160 

0.$16 

400 

165 

0.552 

405 

0.906 

170 

0*652 

410 

0.906 

175 

0.612 

420 

0.906 

180 

in 

• 

© 

440 

0*906 

185 

0.520 

460 

0.906 

190 

0.634 

480 

0.906 

195 

0.662 

500 

0.906 

200 

0.697 

515 

205 

0.527 

520 

210 

0.752 

565 

0*894 

Y/D 

NORTH 


0.606 


0.669 


0.606 


0.656 


0.674 


0.726 


0*787 

0<787 
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MODEL  TEST  34 


BED  MATERIAL— -OTTAWA  SAND 


TEMPERATURE«76.0*  F 


TOTAL  AMPLITUDE*^. 12  IN  D=0.501  IN  o<  *0 

L/D  »  16 


CYCLES  Y/D 


SOUTH 

NORTH 

0 

0.063 

2 

0.227 

5 

0*245 

10 

0.110 

15 

0*263 

20 

0.331 

25 

C.331 

30 

0.276 

95 

0.346 

40 

0.346 

45 

0.346 

50 

0.346 

55 

0.346 

60 

0.346 

65 

0.292 

70 

0.292 

75 

0.339 

80 

0.339 

85 

0.339 

90 

0.339 

95 

0.339 

100 

0.339 

105 

0.292 

110 

0.339 

115 

0*339 

120 

0.339 

125 

0.339 

130 

0.339 

135 

0.339 

CYCLES  Y/D 

SOUTH  NORTH 

140  0.339 

143  0.339 

150  0.345 

135  0.292 

160  0.346 

165  0.346 

170  0.346 

175  0.346 

180  0.346 

185  0.292 

190  0.354 

195  0*354 

200  0*354 

205  0.354 

210  0*354 

215  0.354 

220  0.366 

225  0*366 

230  0.366 

235  0.366 

240  0.366 

245  0.385 

250  0.385 

255  0.385 

260  0.385 

265  0.385 

270  .  0.385 

275  0.385 

310  0.385 

315  0.385 
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MODEL  TEST  37 


BED  MATERIAL — OTTAWA  SAND  TEMPERATURE-75 . 0°  F 

TOTAL  AMPLITUDE *11. 99  IN  D-0.501  IN  04  =0° 


CYCLES  Y/D 


SOUTH 

NORTH 

0 

0.022 

0*022 

2 

0.335 

0.335 

5 

0.367 

10 

0.374 

15 

6.382 

26 

0,524 

25 

0.303 

30 

0.610 

35 

0,287 

40 

0.618 

45 

0.602 

50 

0.397 

55 

0.653 

CYCLES 

Y/D 

SOUTH 

NORTH 

60 

0.378 

65 

0*658 

70 

0*429 

75 

0.607 

80 

0*462 

85 

0*618 

90 

0.488 

95 

0*626 

100 

0.528 

105 

0.650 

110 

0.517 

115 

0.682 

120 

0.532 

125 

0.768 
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MODEL  TEST  38 

8ED  MATERIAL— ‘OTTAWA  SAND  TEMPERATURE=73 . 5°  F 

TOTAL  AMPLITUDE*! 2 • 26  IN  D*0.501  IN  <*0  *0° 

L/D  =  8 (ROUNDED  ENDS) 

CYCLES  Y/D  CYCLES  Y/D 


SOUTH 

NORTH 

SOUTH 

NORTH 

0 

O.064 

170 

0.905 

2 

0.205 

180 

0.913 

5 

0.386 

195 

0.635 

10 

0.414 

200 

0.937 

15 

0*414 

205 

0.953 

20 

0.414 

210 

0.965 

25 

0.414 

215 

0.976 

30 

0.474 

220 

0.630 

35 

0.552 

230 

0.630 

40 

0.442 

235 

0.976 

45 

0.591 

260 

0.607 

50 

0.604 

290 

0.976 

55 

0*488 

293 

0.607 

60 

0.630 

300 

0.984 

65 

0.481 

305 

0.662 

70 

0.685 

315 

0.999 

75 

0.544 

330 

1.000 

80 

0.682 

335 

0*662 

85 

0.701 

355 

0.992 

90 

0.583 

360 

0.653 

95 

0.717 

365 

0.653 

100 

0.717 

370 

0.677 

105 

0.613 

375 

0.945 

110 

0.787 

380 

0.705 

115 

0.803 

395 

0.705 

120 

0.606 

400 

0.945 

125 

0.831 

415 

0.913 

130 

0.843 

420 

0.756 

135 

0.843 

425 

0.910 

140 

0.863 

435 

0.910 

145 

0.614 

440 

0.764 

150 

0.874 

455 

0.897 

155 

0.635 

465 

0.913 

160 

0.997 

470 

0.788 

160 

0.897 

470 

0.788 
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MODEL  TEST  29 

BED  MATERIAL— OTTAWA  SAND  TEMPERATURE® 74. 5°  F 

TOTAL  AMPLITUDE*^. 16  IN  D=0.501  IN  <V  =0° 


CYCLES 

'  Y/D 

SOUTH 

north 

0 

0.245 

0*145 

2 

0.250 

0.250 

5 

0.286 

0*286 

10 

0.354 

15 

0,321 

0.321 

20 

0,335 

0.335 

25 

0.370 

30 

0.4Q2 

0.402 

35 

0.468 

40 

0.419 

45 

0.488 

50 

0,433 

55 

0.535 

60 

0.449 

65 

0.578 

70 

0.457 

75 

0.591 

80 

0,508 

55 

0.559 

0.559 

90 

0.559 

0.559 

95 

0.622 

100 

0,528 

105 

0,646 

110 

0.551 

115 

0.661 

120 

■'  ’  ' 

0.571 

125 

0.646 

130 

0*667 

135 

0.587 

140 

0.701 

145 

0.615 

150 

0.693 

155 

0.693 

CYCLES  Y/D 


SOUTH 

NORTH 

160 

0.725 

0.725 

165 

0«674 

170 

0.740 

175 

0.748 

0.748 

180 

0.748 

0*748 

185 

0.772 

0.772 

190 

0.858 

195 

0.736 

200 

0,882 

205 

0.811 

210 

0.835 

0.835 

215 

0.835 

0.835 

220 

0.851  . 

0.851 

225 

0.886 

230 

0.792 

235 

0.886 

240 

0.891 

245 

0.882 

250 

0.827 

255 

0.966 

260 

0.795 

265 

0.976 

270 

0.799 

275 

0.968 

280 

0.795 

285 

0.968 

290 

0.811 

295 

0.976 

300 

0.811 

305 

0.976 

310 

0.811 

315 

1.000 

320 

0.732 
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MODEL  TEST  41 


BED 

MATERIAL — GLASS 

BEADS 

TEMPERATURE»75  .  tf  F 

TOTAL 

AMPLITUDE=33.80 

IN  D=3 .45 

IN  o <  =0° 

O 

CYCLES 

Y/D 

SOUTH  NORTH 

CYCLES  Y/D 

SOUTH  NORTH 

0 

0*068  0.032 

60  0.792 

5 

0.243 

65  0.801 

10 

0.318 

70  0.811 

15 

0.464 

75  0.844 

20 

0.540 

85  0.909 

25 

0.626 

90  0.932 

30 

0*649 

95  0.959 

35 

0.611 

100  0.989 

40 

0.796 

105  1.012 

45 

0.768 

110  1.023 

50 

0.776 

115  1.041 

55 

0.782 

120  1.064 

126  1.072  1.103 
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MODEL  TEST  42 

BED  MATERIAL — GLASS  BEADS  TEMPERATURE=76» 2*  E 


TOTAL  AMPLITUDE =33 *90  IN  D= 
CYCLES  Y/D 


SOUTH 

NORTH 

0 

0.030 

0.044 

2 

0.182 

5 

0.264 

9.5 

0.372 

15 

0.447 

20 

0.546 

25 

0.550 

30 

0.667 

35 

0.676 

40 

0.692 

45 

0.784 

50 

0.78$ 

•  45  IN  o^o  =0* 

CYCLES  Y/D 

SOUTH  NORTH 

55  0.791 

60  0.792 

65  0.826 

70  0.861 

75  0.886 

80  0.903 

85  0.921 

90  0.942 

95  0.963 

100  0.992 

105  1 .016 

110  1.079  1.107 
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MODEk  TEST  43 

BED  MATERIAL— GLASS  BEADS  TEMPERATURE=74*  5°  F 


TOTAL  AMPLITUDE=33*80  IN  D=3.45  IN  c*  *0° 

© 


CYCLES 

Y/D 

SOUTH 

CYCLES 

Y/D 

SOUTH 

1.92 

0.206 

24.65 

0.582 

3.15 

0.247 

27.40 

0.587 

4.11 

0.284 

30.15 

0.638 

5.03 

0.308 

32.85 

0.639 

6 .04 

0.321 

35.60 

0.653 

7.21 

0.326 

38.35 

0.710 

8.04 

0.326 

41,15 

0.762 

9.28 

0.352 

44.10 

0.762 

11.77 

0.420 

46.60 

0.765 

11.95 

0.441 

49.35 

0.768 

13.21 

0.452 

52.20 

0.771 

14.25 

0.456 

55.20 

0.778 

15.21 

0.459 

59.30 

0.793 

16.30 

0.465 

64.30 

0.799 

17,62 

0.512 

68.90 

0.820 

19.00 

0.560 

74.80 

0.862 

19.78 

0.562 

79.60 

0.888 

20.82 

0.564 

86.50 

0.927 

22.60 

0.572 

93.40 

0.950 

95.80 

0.963 
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MODEL  TEST  44 


BED  MATERIAL — GLASS  BEADS 
TOTAL  AMPLITUDE=33.50  IN  D= 
CYCLES  Y/D 


SOUTH 

NORTH 

0 

0.195 

0.204 

2 

0.195 

5 

0.240 

10 

0.380 

15 

0.447 

20 

0.525 

25 

0.560 

30 

0.615 

35 

0.664 

40 

0*695 

TEMPERATURE=75  •  0*  F 

.45  IN  0^*0° 

CYCLES  Y/D 

SOUTH  NORTH 

45  0.733 

50  0.775 

55  0.813 

60  0.884 

65  0.876 

70  0.905 

75  0.923 

80  0.942 

85  0.959 

90  0.979 

95  0.993 
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MODEL  TEST  45 


BED  MATERIAL-GLASS  BEADS 
TOTAL  AMPLITUDES. 60  IN  D« 
CYCLES  Y/D 


SOUTH 

NORTH 

0 

0.287 

0.250 

6 

0.303 

11 

0.446 

16 

0.492 

21 

0.536 

26 

0.618 

31 

0.653 

36 

0.679 

41 

0.707 

46 

0.748 

TEMPERATURE *7 5 . 8  F 

•  45  IN  =0°  ' 

CYCLES  Y/D 

SOUTH  NORTH 

51  0.788 

56  0.810 

61  0.856 

66  0.893 

71  0.914 

76  0.926 

81  0.943 

86  0.956 

92  0.999 

96  1*011 

101  1.023 
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MODEL  TEST  46 


BED  MATERIAL-GLASS  BEADS 

TOTAL  AMPLITUDE=33«40  IN  D* 

CYCLES  Y/D 

SOUTH  NORTH 

0.526  0.482 

0.354 
0.5  54 
0.554 
0.554 
0.554 
0.554 
0.662 
0.676 
0.708 
0.742 
0.776 


TEMPERATURE=75 . 5  F 

.45  IN  »0° 

CYCLES  Y/D 

SOUTH  NORTH 

55  0.792 

60  0.805 

65  0.829 

70  0.850 

75  0.866 

80  0*886 

85  0.904 

90  0.914 

95  0.925 

100  0.936 

105  0.951 

110 
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MODEL  TEST  47 


BED  MATERIAL — GLASS  BEADS  TEMPERATURE=76 . 8°  F 


TOTAL 

AMPLITUDE *33 .60  IN 

D= 1 • 702  IN  c* 

o 

O 

=0 

CYCLES 

Y/D 

CYCLES 

Y/D 

SOUTH  NORTH 

SOUTH  NORTH 

0 

0.217  0.205 

16 

0.739 

4 

0.371 

20 

0.787 

8 

0.554 

24 

0.905 

12 

0.625 

28 

0*954 

32 

0.984 
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MODEL  TEST  50 

BED  MATERIAL— GLASS  BEADS  TEMPERATURES. 6° F 

TOTAL  AMPLITUDES. 10  IN  D=1.702  IN  K,  =0' 


CYCLES 

Y/D 

CYCLES 

Y/D 

SOUTH 

SOUTH 

0 

0.058 

I  2 

0.743 

4 

0.460 

16 

0.801 

8 

0.616 

20 

0.861 

24 

0.97  7 
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DISCUSSION  OF  RESULTS 

The  model-test  program  was  executed  at  two  values  of  the  sediment 
Froude  number,  U  / >/(s-l)  gd;  namely,  2.8  for  Tests  25-39,  inclusive  and 
11.2  for  Tests  5~l8>  4l-47j  and  50-  Model  tests  were  impossible  to  per¬ 
form  at  intermediate  values  of  the  sediment  Froude  number,  F,  because  of 
the  formation  of  ripples  on  the  bed.  These  bed  ripples  are  of  the  same 
order  of  size  as  the  models  and,  in  fact,  the  cylindrical  models  become 
part  of  the  ripple  system  with  the  axis  of  the  cylinder  coinciding  with 
a  ripple  crest.  This  condition  is  illustrated  in  Figure  4  for  which  the 
value  of  F  was  5-9-  When  the  model  is  part  of  the  ripple  system  no  scour 
occurs.  In  contrast,  when  the  bed  ripples  are  nonexistent  or  insignificant 
in  height  a  scour  hole  develops  around  the  model  as  shown  in  Figure  5  for 
which  the  value  of  F  was  11.2. 

The  lower  value  of  F,  that  is  2.8,  was  selected  as  the  maximum  value 
at  which  model  tests  could  be  made  without  the  formation  of  ripples  on 
the  bed.  This  value  was  selected  by  considering  the  critical  scour  para¬ 
meter  (Appendix  A) 

C  =  =  -  (2) 

v/(s-l)gd(tan  4>  cos'*  +  sin«*  J” 

in  which  ot  is  the  angle  of  inclination  of  the  bed  with  the  horizontal. 

The  value  of  C  for  incipient  motion  of  bed-particles  in  oscillatory  flow 
is  about  3-5  (Table  3,  Appendix  A).  The  angle  of  repose  of  the  Ottawa 
sand  used  in  Tests  25-39)  inclusive,  is  32-5  degrees.  For  an  originally 
level  bed  <=<  is  zero.  Substituting  these  values  into  equation  (2)  the 
value  of  F  of  2.8  is  obtained. 


CONFIDENTIAL 


CONFIDENTIAL 

50 


Figure  4.  Model  Incorporated  into  Ripple  System  (F  =  5.9). 
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The  criterion  that  F  be  equal  to  or  less  than  2.8  was  found  to  be 


inadequate.  In  the  following  table  are  listed  the  model  tests  and  the 
pertinent  characteristics  which  were  performed  for  the  purpose  of 
determining  the  conditions  of  ripple  formation. 


Test  No. 

D(in) 

d(mm) 

D/d 

u 

m 

y(s-i)gd 

Bed  Condition 
Around  Model 

1 

1.002 

O.297 

85.8 

1-39 

scoured 

2 

1.002 

0.297 

85.8 

2.85 

rippled  from  beginning 

3 

1.702 

0.297 

145.6 

2.12 

rippled  @  1340  cycles 

20 

1.002 

0.585 

43.6 

2.l6 

rippled  @  300  cycles 

21 

1.002 

0.585 

43.6 

2.72 

rippled  @  150  cycles 

22 

1.002 

O.585 

43.6 

3-21 

rippled  @  120  cycles 

23 

1.702 

O.585 

74.0 

3-15 

rippled  from  beginning 

24 

0.501 

0.585 

21.8 

2.78 

scoured 

48 

0.501 

0.297 

42.9 

3-25 

rippled  @  35  cycles 

The  tests  results  shown 

above  indicate  that 

the  bed  is 

quasi-stable  in  that 

bed  ripples 

will  form  at  values  of 

F  less  than  for  the 

incipient-motion 

condition; 

for  example, 

in  Test  20. 

Ripple 

inception  is  definitely  a 

function  of 

F  as  shown  by  Tests  20, 

21,  and 

22  in  which  the  time  for  ripple 

inception  to  occur  decreased  as  the  value  of  F  was  increased.  In  addition, 
the  magnitude  of  the  disturbance  is  quite  significant  as  evidenced  by 
Tests  2  and  24  for  which  the  values  of  F  were  nearly  equal.  Ripples  existed 
from  the  beginning  of  Test  2  but  did  not  exist  in  Test  24.  The  principal 
difference  between  these  two  runs  is  in  the  magnitude  of  the  disturbance 
as  evidenced  by  the  cylinder  diameter-to-partiele  diameter  ratio.  Mo 
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further  investigation  of  the  role  of  disturbance  magnitude  in  ripple 
inception  was  made.  The  decision  was  made  to  conduct  some  of  the  model 
tests  with  the  0.501-in  diameter  model  on  the  Ottawa-sand  bed  at  a 
value  of  F  of  2.8. 

The  higher  value  of  F,  that  is  11.2,  was  selected  as  the  minimum 
value  at  which  the  bed  ripples  are  sheared  off.  In  this  range  the 
uppermost  particles  of  the  bed  are  in  motion  over  the  entire  bed  and 
ripples  disappear.  Inman  (2)  presents  results  from  the  studies  of 
Manohar  (3)  and  Menard  (4)  showing  that  the  velocity  at  which  ripples 
disappear  increased  with  increasing  grain  size.  This  observation  was 
verified  in  the  present  experiments  in  that  ripples  would  disappear  or 
be  insignificant  at  a  maximum  velocity,  U^,  of  2.4  fps  with  the  0.297 -mm 
diameter  bed  material  but  the  ripples  would  not  disappear  at  a  value  of 
U  of  2.7  fps  with  the  O.585  nan  diameter  bed  material.  Since  U  of  2.7 
fps  was  the  maximum  attainable  in  the  U-tube,  the  tests  in  which  the 
ripples  disappeared  were  limited  to  the  smaller  bed  material  with  a 
value  of  F  of  11.2. 

Phenomenological  Discussion  of  Tests  with  F  of  2.8 

The  following  discussion  is  presented  chronologically.  In  all  tests, 
the  starting  transients  were  eliminated  by  utilizing  an  initial  blowdown 
in  the  East  vertical  leg  of  the  U-tube. 

Immediately  upon  the  start  of  a  run  scour  holes  begin  to  form  at  the 
ends  of  the  cylinder.  These  scour  holes  deepen  and  expand  with  much  of 
the  scoured  material  being  carried  toward  the  middle  of  the  cylinder 
forming  a  central  ridge  which  decreases  in  height  away  from  the  model. 
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When  the  scour-hole  bottom  is  below  the  bottom  of  the  model,  material 
slides  out  from  under  the  model  forming  a  flow  passage  under  the 
cylinder.  Once  the  flow  passage  under  the  model  is  started,  material 
is  carried  under  the  cylinder  and  these  flow  passages  increase  in  size 
by  lateral  movement  toward  the  center.  The  central  support  upon  which 
the  cylinder  rests  diminishes  in  size  until  the  support  is  inadequate 
at  which  time  the  cylinder  rocks  on  the  support  and  settles  rapidly. 

Upon  settling  the  supporting  base  is  broadened  only  to  be  narrowed 
again  by  scour  under  the  cylinder.  The  burial  proceeds  in  this  manner 
settling  spasmodically  with  ever  increasing  time  intervals  between 
settlement  increments .  The  model  settles  in  an  essentially  level  manner 
upon  collapse  of  the  central  support.  The  largest  observed  deviation 
from  the  horizontal  was  2 6  degrees. 

The  material  which  is  removed  from  under  the  cylinder  must  obviously 
be  disposed  of  in  some  manner.  The  large  eddies  formed  at  the  ends  of 
the  cylinder  appear  to  be  the  main  mechanism  for  excavating  and  removing 
material  scoured  from  under  the  cylinder.  In  the  tests  with  F  of  2.8,  the 
material  was  moved  by  rolling  along  the  bed.  This  type  of  movement  is 
called  bed-load  movement.  The  fully  developed  scour  hole  resembles  an 
inverted  frustum  with  the  peripheral  surface  being  inclined  at  the  angle 
of  repose  of  the  bed  material. 

In  summary,  the  settlement  can  be  described  as  a  destruction  of 
the  cylinder  support  by  scour  under  the  cylinder.  Simultaneously  the 
eddy  pattern  created  as  a  result  of  flow  separation  is  the  mechanism 
for  the  enlargement  of  the  surrounding  scour  hole.  Both  scour  in  the 
vicinity  of  the  model  and  under  the  model  are  necessary  for  settlement 
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to  occur.  In  fact  the  rate  at  which  material  can  he  scoured  from  under 
the  cylinder  is  limited  by  the  rate  at  which  material  can  he  carried 
out  of  the  surrounding  scour  hole. 

Phenomenological  Discussion  of  Tests  with  F  of  11.2 

As  before  the  discussion  is  chronological. 

Upon  starting,  the  model  rolls  to  and  fro  in  the  first  cycle  and 
settles  rapidly.  No  rolling  occurs  after  the  first  cycle.  Visual 
observation  of  the  progress  of  scour  is  difficult  since  the  removed 
material  leaves  the  scour  hole  in  suspension  as  shown  in  Figure  6, 

There  is  no  reason  to  believe  that  the  basic  process  of  settlement  is 
any  different  in  the  tests  in  which  F  is  11.2  than  those  in  which  F  is 
2.8.  The  difference  is  solely  in  the  mechanism  of  removal;  that  is, 
as  suspended  load  in  the  former  and  as  bed  load  in  the  latter.  No 
difference  could  be  detected  in  the  scour-hole  geometry  between  the 
two  sets  of  tests. 

The  cylinder  rests  on  a  fulcrum  at  the  instant  of  settling  which 
permits  the  cylinder  not  only  to  rock  but  also  to  turn.  When  the  initial 
orientation  q  is  different  from  zero  the  cylinder  also  pivots  during 
settlement.  The  cylinder  continues  to  pivot  in  increments  until  the 
axis  is  parallel  to  the  wave  crest. 
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Figure  6.  Removal  by  Suspended  Transport  (F  =  11.2) 
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ANALYSIS  OF  RESULTS 

The  model  study  was  severely  limited  by  the  existence  of  ripples 
of  the  same  order  in  size  as  the  model.  Only  two  values  of  the  sediment 
Froude  number,  F,  could  be  obtained.  The  higher  value,  F  =  11.2,  cor¬ 
responds  to  a  gravity  wave  system  generated  by  a  severe  storm j  whereas, 
the  lower  value,  F  =  2.8,  corresponds  to  a  wave  system  generated  in  a 
moderate  sea.  Since  the  experimental  method  could  not  be  employed  for 
the  intermediate  range  of  F,  the  writers  were  forced  to  make  a  detailed 
analysis  of  the  results  of  other  scour  studies.  The  results  of  two 
excellent  studies,  Laursen  (5  }  and  Rouse  (6),  of  jet  scour  were  analyzed 
(Appendix  A) .  In  neither  of  these  studies  was  ripple  formation  a  factor. 
Consequently  both  Laursen  and  Rouse  were  able  to  vary  systematically  the 
sediment  Froude  number.  The  sediment-transport  functions  (Appendix  A, 
equations  11  and  l6)  form  the  basis  for  interpolating  and  extrapolating 
these  model  test  results  throughout  the  entire  region  of  interest. 
Scour-Hole  Geometry 

The  scour  hole  formed  around  the  cylinder  is  three-dimensional  in 
form  as  shown  in  Figure  5*  After  Test  47,  the  scour  hole  was  mapped  by 
using  a  vertical  point  gage  to  determine  the  elevation.  The  point  gage 
was  located  in  various  positions  in  the  horizontal  plane.  From  these 
measurements  a  contour  map  of  the  scour  hole  was  prepared  as  shown  in 
Figure  7 •  The  scour-hole  can  be  represented  as  an  inverted  frustum  with 
the  side  slope  equal  to  the  angle  of  repose,  <p  ,  of  the  bed  material. 
Using  the  volume  computed  from  Figure  7,  the  diameter  of  the  bottom  sur¬ 
face  (bottom  of  the  scour  hole)  of  the  frustum  is  L  +  0.24  D  in  which 

CONFIDENTIAL 


57 


CONFIDENTIAL 


i 


CONFIDENTIAL 


L  and  D  are  the  cylinder  length  and  diameter,  respectively.  Using 
this  representation  the  volume  of  the  scour  hole  ¥■,  is  related  to 
the  depth  as  follows. 


Since  L/D  is  constant  at  a  value  of  4,  equation  (3)  can  he  simplified  to 


The  rate  that  bed  material  is  removed  from  the  scour  hole  is  simply 
the  time  rate  of  change  of  ¥•  with  respect  to  time,  t,  that  is 


n  -  ^ 

Qso  dt  (5) 

in  which  is  the  sediment  discharge  out  of  the  scour  hole. 

Performing  the  operation  indicated  in  equation  (5)  on  equation  (4). 


Making  equation  (6)  dimensionless. 


Qso  (  3-144  13-4  l4.l4  \  d(  ys/D) 

Um  yg2  \tan2c()  (  yg/D')tan4>  ( jjv  )2  J  d^t/D)  (7) 
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All  of  the  quantities  on  the  RHS  of  equation  (7)  can  be  evaluated  from 

experimental  results.  As  an  aid  in  determining  the  slope,  d(yc  /d) / 

d ( U^t/D ) ,  the  data  v/ere  plotted  on  logarithmic  paper.  Only  data  for 

which  the  scour  hole  was  well-formed,  that  is  yg/D  >0-5,  were  used 

2 

in  computing  the  transport  function,  Q_/u  yg 

As  explained  previously  the  model  tests  were  severely  limited  in 
regard  to  the  variation  of  the  sediment  Froude  number,  U  /  \/(s-l)  gd  . 
This  limitation  was  overcome  by  analyzing  the  scour  experiments  of 
Laursen  (5),  Rouse  (6),  and  Ahmad  (7)  (Appendix  A).  From  these 
analyses,  the  writers  concluded  that  the  sediment  transport  functions 
of  scour  are  of  the  following  form 


for  suspended  transport  and 

QSO  _  8  /  \  m 

F  L  ) 

m  s 

for  bed-load  transport.  The  exponents  of  the  F  terms  were  determined  from 
the  experiments  of  others.  The  exponents  of  yg/D  were  determined  using 
equation  (7)  and  model-test  data. 

Tiie  resulting  transport  functions  for  a  cylinder  lying  on  the  sea 
bed  under  an  oscillatory  wave  are  shown  graphically  in  Figure  8  for 
suspended  transport  and  in  Figure  9  for  bed-load  transport.  These 
functions  are 
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Figure  8-  Suspended-Load  Transport  Function  Tor  a  Horizontal 
Cylinder  ( L/D  =  4 ) . 


CONFIDENTIAL 

6l 


CONFIDENTIAL 


-3,3  I _ _ _ l_^ _ _ _ _ I _ I _ I - - 1 - 1 - 1 - J 

15  3.4  3.7  38  39  4.0  4.1  4.2  4.3  4.4 


log 


8 


F 


-3 


Figure  9.  Bed-Load  Transport  Function  for  a  Horizontal 
Cylinder  (L/D  =  4). 


CONFIDENTIAL 

62 


CONFIDENTIAL 


1/2 


u  y 
m  s 


O.Ol 


< )  ^ 


-3 


for  suspended  transport  and 


(8) 


Q 


SO 


=  1.2  (lO'7)  F8  ^  — ^ 


(9) 


for  bed-load  transport.  Classification  as  to  the  mode  of  transport  is 
required  since  the  transport  was  by  bed-load  movement  during  model  tests 
with  F  =  2.8  (Tests  25-39,  inclusive)  and  the  transport  was  by  suspended 
load  movement  with  F  =  11.2  (Tests  5-13,  inclusive,  4l-47,  inclusive, 
and  50 ) . 

Settlement  Functions 

Having  obtained  the  sediment-transport  functions,  equations  (8) 
and  (9) ,  equation  (6)  can  be  integrated  to  obtain  the  settlement  functions. 
The  reason  for  analyzing  the  problem  in  this  apparently  roundabout  manner 
is  that  prior  history  is  removed  by  evaluating  a  rate  function.  Thus  the 
constant  of  integration  obtained  from  the  settlement-differential  equation, 
equation  (6),  will  be  a  function  of  the  initial  conditions.  Substituting 
equation  (8)  into  equation  (6)  and  integrating 


(io'2)f 


+  K 


0.786  , 

1  i  1 

fh) 

J 

1  +  7.07 | 

fl±  \ 

,  2  1 
tan  <p 

4  D  / 

1  +  1 
tan 

l  D  / 

U  /  (10) 

in  which  Kg  is  the  constant  of  integration  to  be  obtained  using  model- 
test  data.  A  similar  expression  is  obtained  upon  substitution  of 
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equation  (9)  into  equation  (6)  and  integrating,  that  is, 


1.2  (10"7)F8 


k 

0,786— 

tan  '  D  / 


Again  the  constant  of  integration,  K^,  can  be  obtained  by  using  the 
model-test  data. 

Initial  Conditions 

The  value  of  the  constant  of  integration,  K  in  equation  (10)  and 

s 

in  equation  (ll),  is  dependent  upon  initial  conditions.  The  obvious 
initial  conditions  are  the  initial  burial,  y  /D,  and  the  initial 
orientation,  oc  ^  .  Two  other  conditions,  initial  development  of  scour 
hole  and  initial  scour  under  the  cylinder,  are  not  initial  conditions 
but  can  be  treated  as  initial  conditions .  The  constants  of  integration 
were  computed  by  substituting  values  of  the  experimental  parameters,  F, 

,  and  d/D  and  of  the  experimentally  determined  parameters,  U^t /l) 
and  yo/D,  into  equations  (10)  and  (ll)  for  each  test.  Only  data  for 
which  the  scour  hole  was  fully  developed,  that  is  y^/D  >  0.5,  were 
used  in  the  determination  of  K.  An  arithmetic  mean  of  the  values  of 
the  K's  for  each  datum  point  was  used  as  the  correct  value  of  K  for 
that  test. 

Initial  Burial  -  The  constant,  K^,  of  equation  (ll)  was  determined 
using  the  method  described  above  using  the  data  Tests  25,  27,  28,  29, 

30,  and  39  for  which  ^  was  zero.  The  value  of  was  found  to  vary 
systematically  with  initial  burial,  y^/D .  The  values  of  obtained 
from  these  six  tests  are  shown  in  Figure  10.  The  constant  of  integration 
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K^,  in  the  settlement  equation,  equation  (ll),  can  be  approximated 


^  =  5-2 


0.63 


(12) 


A  similar  procedure  was  applied  to  the  data  of  Tests  8,  13, 

4l-45,  47,  and  50  in  order  to  determine  the  constant  of  integration, 

Ks,  in  the  settlement  equation,  equation  (10),  when  the  material 
was  removed  as  suspended  load .  No  variation  in  Kg  could  be  detected 
with  values  of  initial  burial,  y  /D,  which  ranged  from  0.037  for 
Test  4l  to  0.269  for  Test  45-  In  fact  for  this  series  of  tests  the 
value  of  K  was  negligible.  Hence  the  value  of  Kg  in  equation  (10) 
was  taken  as  zero  for  <Xq  =  0  irrespective  of  initial  burial. 

Initial  Orientation  -  Whenever  the  angle  <x.  between  the  cylinder 
axis  is  different  from  0  or  90  degrees,  the  flow  exerts  a  torque  on 
the  cylinder  tending  to  decrease  the  angle .  A  cylinder  placed  at  an 
initial  orientation  of  Q  pivots  in  increments  until  oc  is  zero. 

The  settlement  of  the  cylinder  occurs  in  increments  as  the  central 
support  under  the  cylinder  is  scoured  away.  The  central  support  under 
the  cylinder  also  serves  as  a  pivot  for  turning.  The  history  of 
turning  is  shown  in  Figure  11  in  which  data  from  Tests  5-7,  10-12,  and 
l4-l8  are  included.  If  oc  were  less  than  60  degrees,  (a)  the  cylinder 
turned  rapidly;  (d)  an  ordinary  scour  hole  developed;  and  (c)  settle¬ 
ment  occurred  at  essentially  the  same  rate  as  a  cylinder  placed  with 

oc  =0.  Therefore,  for  values  of  OC  less  than  60  degrees,  the 
0  0 

constants  of  integration  required  no  adjustment  for  initial  orientation. 
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Figure  11.  History  of  Turning. 
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With  values  of  o£q  greater  than  60  degrees  an  abnormal  scour  hole 
developed  for  which  the  analysis  based  upon  geometrically  similar  scour 
patterns  would  not  apply.  As  a  consequence,  the  analysis  and  the  results 
presented  in  this  report  apply  only  if  is  less  than  60  degrees.  This 

limitation  is  not  serious  because  2/3  of  the  cases  of  random  placement 
can  be  analyzed  and  because  the  settlement  rate  is  not  drastically  dif¬ 
ferent  for  the  remaining  l/3  of  the  cases. 

Initial  Development  of  the  Scour  Hole  -  The  scour  hole  begins  as 

two  isolated  scour  holes  at  the  ends  of  the  cylinder.  The  scour  hole 

continues  to  change  in  geometric  configuration  until  the  cylinder  has 

settled  about  one-half  a  diameter,  that  is,  until  y  /D  =  0.5.  The 

s 

scour  hole  is  geometrically  similar  with  further  settlement  having  the 
geometric  features  shown  in  Figures  5  and  7-  The  constant  of  integration, 
Kg  in  equation  (10)  and  in  equation  (ll)  is  computed  from  experimental 
measurements  taken  after  the  scour  hole  is  well -developed .  Thus  the 
value  of  the  constant  of  integration  determined  in  this  manner  reflects 
the  effect  of  the  development  stage. 

Initial  Scour  Under  the  Cylinder  -  Settlement  of  the  cylinder  occurs 
in  steps .  Scour  occurs  under  the  cylinder  decreasing  the  width  of  the 
central  support  until  the  central  support  collapses .  Then  the  process 
is  repeated.  Since  scour  under  the  cylinder  is  the  direct  cause  of 
settlement,  an  exploratory  study  was  conducted  of  flow  under  a  partially 
imbedded  cylinder.  The  results  are  presented  in  Appendix  B.  At  a  later 
date,  the  conclusion  was  reached  that  the  settlement  rate  was  limited 
to  the  rate  at  which  material  could  be  removed  from  the  scour  hole  rather 
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than  from  under  the  cylinder.  Consequently ,  the  final  analysis  presented 

in  this  report  is  based  entirely  upon  the  scour-hole  geometry  and  rate  of 

removal  of  material  from  the  scour  hole. 

Summary  -  Two  initial  conditions  have  been  analyzed  in  respect 

to  the  effect  of  these  conditions  upon  the  value  of  the  constant  of 

integration  in  the  settlement  equation.  The  initial  burial,  y^/l), 

was  found  to  influence  the  constant  of  integration  when  the  bed 

material  was  removed  by  bed-load  transport.  On  the  other  hand  no 

effect  of  initial  burial  could  be  discerned  when  the  bed  material  was 

removed  by  suspended- load  transport.  The  constant  of  integration  was 

not  influenced  by  the  initial  orientation,  «  ,  if  <=^  was  less  than 

o  o 

60  degrees .  Even  for  values  of  c*  ^  greater  than  60  degrees  following 
settlement  functions  are  tolerable. 

Inserting  the  constant  of  integration  into  equation  (10), 


0.786 

2 

tan 


is 

tanijX  D 


(10a) 


Equation  (10a)  is  the  settlement  function  for  a  horizontal  cylinder 
of  L/D  =  4  when  the  removal  of  bed  material  is  by  suspended-load 
transport.  All  of  the  experimentally  determined  data  from  Tests  5-18, 
4l-47,  and  50  are  shown  in  Figure  12.  Equation  (10a)  is  also  shown  in 
Figure  12  in  which  a  value  of  4>  of  24  degrees  was  used.  On  account 
of  the  difficulty  of  solving  equation  (lOa)  for  a  value  of  y  /D, 
Figure  13  has  been  prepared  as  a  graphical  aid  for  problem  solving. 
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Figure  12.  Experimental  Results  for  Suspended-Load.  Transport  Tests. 
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Figure  13.  Settlement  Function. 
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Inserting  the  constant  of  integration,  equation  (12)  into 
equation  (ll), 


L.2  (10'7)F8(^  — j  +3.2^—^ 


°'63  0.786  /  ys  \\  ^5  y  3 

o  /  tan  4>^  D  ' 


tan  41  V  D 


7.07 


& 


(lla) 


Equation  (lla)  is  the  settlement  function  for  a  horizontal  cylinder  of 
L/d  =  4  when  the  removal  of  bed  material  is  by  bed-load  transport.  All 
of  the  experimentally  determined  data  from  Tests  25 -3°,  37,  and  39  are 
shown  in  Figure  l4.  Equation  (lla)  is  also  shown  in  Figure  l4  in  which 
a  value  of  (f)  of  32.5  degrees  was  used.  On  account  of  the  difficulty  of 
solving  equation  (lla)  for  a  value  of  yg/D,  Figure  13  has  been  prepared 
as  a  graphical  aid  for  problem  solving. 

Range  of  Applicability 

The  question  arises  as  to  which  of  the  settlement  functions,  equa¬ 
tion  (lla),  is  applicable  in  a  given  situation.  The  limit  between  bed 
load  transport  is  found  by  equating  the  RHS  of  equation  (8)  to  the  REIS 
of  equation  (9)  and  then  solving  for  the  limiting  value  of  F.  Per¬ 
forming  these  operations 


F( limit) 


5.04  (d/D)1/1^ 


(13) 


The  limiting  value  of 
sediment  diameter,  d, 
table. 


F>  U  /,/(s-l)gd  ,  is  weakly  dependent  upon  the 
/ 

and  the  model  diameter,  D,  as  shown  in  the  following 
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Figure  l4. 


Experimental  Results  for  Bed-Load.  Transport. 
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d(mm) 

D(in) 

F(limit) 

0.2 

18 

2.9 

0.4 

18 

3-0 

0.6 

18 

3.1 

If  the  value  of  F  is  greater  than  the  value  in  the  third  column, 
suspended-load  transport  predominates  and  the  settlement  equation, 
equation  (lla),  applies.  Conversely  if  the  value  of  F  is  less  than 
the  value  in  the  third  column  of  the  table,  bed-load  transport  pre¬ 
dominated  and  the  settlement  equation,  equation  (10a),  applies. 

As  discussed  previously,  neither  settlement  equation,  equation  (lOa) 
nor  equation  (lla),  applies  when  the  cylinder  becomes  part  of  the  ripple 
system.  In  fact  settlement  ceases  when  the  cylinder  becomes  part  of  the 
ripple  system.  The  question  then  arises  as  to  whether  and  under  what 
conditions  would  a  prototype  mine  become  part  of  the  ripple  system.  A 
rational  criterion  can  be  formulated  based  upon  the  concept  that  if  the 
ratio  of  the  ripple  amplitude,  q  ,  to  the  cylinder  diameter,  D,  exceeds 
a  certain  value,  II,  then  the  cylinder  will  be  a  part  of  the  ripple 
system,  that  is  if  ^  /D  K.  From  Inman  (2) 

Tv/  6 

In  which  Tv.  is  the  wave  length  of  the  ripples .  Also  from  Inman '  s  results 
the  following  empirical  equation  can  be  derived 

"X  s=  1200  d 
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Combining  equations  (l4)  and  ( 15 )  with  the  definition  of  H 
N  =  £  -  200  (  |  ) 

A  reasonable  estimate  is  that  if  S  4  l/4  bed  ripples  will  not 
influence  the  settlement.  Using  this  criterion  in  equation  (l6)  for 
an  l8-in  diameter  cylinder,  no  influence  of  bed  ripples  is  anticipated 
if  the  sediment  diameter  is  less  than  0.57  mm  and  if  there  is  no  initial 
burial. 

If  the  mine  has  been  subjected  to  a  prior  storm  in  which  settlement 
occurred  and  after  which  the  scour  hole  filled,  the  profusion  above  the 
bed  would  be  less  than  D  as  assumed  in  the  above  example.  If  the  mine 
had  previously  buried  to  ys/D  =  0.5,  the  profusion  would  te  9  in  rather 
than  18  in.  Again  one  might  anticipate  that,  if  the  ripple  height  were 
one  fourth  of  the  protruded  height,  the  mine  would  become  part  of  the 
ripple  system.  Since  equations  (l4),  (15),  and  (l6)  are  all  linear  re¬ 
lations  bed- material  diameter  would  simply  be  one  half  of  0.57  mm  or 
0.28  mm.  In  other  words,  if  a  mine  had  buried  halfway  as  the  result  of 
a  previous  storm,  if  the  scour  hole  were  filled,  if  the  bed- material 
diameter  were  greater  than  0.28  mm,  and  if  F  were  in  the  range  of  2  to  10 
the  mine  could  be  expected  to  become  part  of  the  ripple  system  with  no 
settlement  during  the  second  storm.  The  implications  are  that  if  a  mine 
does  not  become  completely  buried  during  the  first  storm  after  placement 
the  mine  is  more  likely  to  become  a  part  of  the  ripple  system  in  suc¬ 
ceeding  storms . 
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The  range  through  which  the  settlement  functions,  equations  (lOa) 
and  (lla)  apply  are  summarized  in  Figure  15 •  The  boundaries  of  the  "no 
settlement1  region  were  not  established  in  this  model  study.  The  upper 
boundary  was  established  from  Manohar's  (3)  experimental  results  for  the 
disappearance  of  ripples.  The  vertical  boundary  at  d  =  0-57  rani  or 
0.28  mm  was  established  as  explained  above.  The  lower  boundary  is  very 
approximately  located  on  the  basis  of  the  incomplete  tests  mentioned  in 
DISCUSS IOW  OF  RESULTS.  The  "no  settlement"  zone  as  shown  in  Figure  15 
is  qualitatively  verified  by  the  observations  reported  by  Salsman  and 
Talbert  (8).  They  report  a  case  in  which  two  mines  were  located  in 
proximity.  The  bed  material  surrounding  one  mine  had  a  median  diameter 
of  0.550  mm  and  the  bed  material  surrounding  the  other  had  a  mean  dia¬ 
meter  of  0.112  mm.  Both  mines  were  subjected  to  the  same  wave  action. 

The  mine  resting  on  the  coarser  bed  material  became  part  of  the  ripple 
system  whereas  a  normal  scour  pattern  developed  around  the  mine  placed 
on  the  finer  material. 

Comparison  with  Prototype 

The  field  tests  conducted  by  the  Navy  Mine  Defense  Laboratory  (USNMDL) 
during  the  winter  of  1962-63  will  be  analyzed  in  order  to  compare  observed 
field  results  with  predictions  based  upon  the  model  tests.  The  mines  were 
dropped  onto  the  sea  bed  and  their  settlement  measured  at  various  intervals. 
Some  mines  were  repositioned  by  divers  to  the  prone  position  (y  /D  =  0) 
after  the  measurements  were  recorded  whereas  other  mines  were  left  in  their 
scour  holes.  Hence  the  effect  of  initial  burial  on  settlement  could  be 
determined.  The  pressure  signature  on  the  sea  bottom  was  recorded  in  order 


CONFIDENTIAL 

76 


CONFIDENTIAL 

SETTLEMENT  (Suspended-load  transport, 
equation  (lOa)  ) 

IV 

1  \ 
i 

1 

lysl/D  =  o.5 

1 

ysl/D  =  0 

NO  SETTLE 

®iNT 

1 
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ripple  system) 
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Equation  (13) 

- - - - 

- -  - 

SETTLEMENT  (Bed -load  transport 

,  equation  (lla)  ) 
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Figure  15.  Range  of  Applicability  of  Settlement  Functions 
for  D  =  l8  in. 
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that  the  wave  amplitude  and  bottom  velocity  could  he  determined.  A 
total  of  three  mines  were  placed  on  the  sea  bed  at  two  sites,  having 
water  depths  of  40  ft  and  60  ft. 

Bed  Materials  -  The  USNMDL  took  core  samples  of  the  bed  materials 
at  the  two  sites  and  ran  size  distribution  analyses.  The  average  size 
of  the  sediment  at  the  bed  is  0.221  mm  for  the  4o-ft  depth  and  O.360  mm 
for  the  60-ft  depth.  The  sand  will  be  assumed  to  have  a  density  cor¬ 
responding  to  that  of  quartz,  i.e.,  in  sea  water,  s  -  2.6o.  The  angle 
of  repose  will  be  assumed  to  be  35°  f°r  the  larger  sand  and  30°  for  the 
smaller. 

Mines  -  A  MK-39  mine  (L  =  88  in,  D  =  22.5  in)  was  dropped  at  the 
40-ft  depth  site  on  21  December  1962,  and  removed  on  27  March  19&3 • 

During  this  period  the  settlement  of  the  mine  was  measured  at  intervals. 

Two  MK-36  mines  (L  =  71  in,  D  =  18.5  in),  referred  to  as  A  and  B  by  USNMDL, 
were  dropped  at  the  60-ft  depth  site  on  12  October  1962 .  Mine  B  was 
repositioned  during  some  of  the  visits  made  for  settlement  measurements 
while  mine  A  was  allowed  to  settle  continually  from  21  December.  The  mine 
at  the  40-ft  depth  will  be  referred  to  as  mine  C  for  convenience.  The 
dimensions  of  the  mines  and  sediment  characteristics  at  the  mine  sites 
are  tabulated  below: 


Mine  Mine  Water  Mean  Sediment  Angle  of 


Mine 

Length  (L) 

Diameter  (D) 

Depth 

Diameter ,  d 

Repose,  <{>  (Assumed) 

A 

71  in 

18.5  in 

60  ft 

O.360  mm 

35  degrees 

B 

71  in 

18.5  in 

6o  ft 

O.360  mm 

35  degrees 

C 

88  in 

22.5  in 

4o  ft 

0.221  mm 

30  degrees 
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Sea  Conditions  -  Personnel  of  the  USNMDL  have  resolved  a  pressure 
spectrum  and  a  bottom  velocity  spectrum  from  recorded  bottom  pressures 
at  the  60-ft  depth  mine  site.  The  velocity  spectrum  data  were  further 
analyzed  and  the  workable  results  along  with  a  letter  of  description 
were  communicated  by  the  USNMDL.  In  the  letter  of  description: 

" .  the  velocity  estimates  represent  the  amplitude  A  in 

Acos  oa  t  (where  to  is  the  frequency  of  the  peak  in  the  spectrum) . 
This  amounts  to  replacing  the  original  random  spectral  process 
by  a  single  frequency  constant  amplitude  process  which  has  the 
same  energy  in  the  root  mean  square  sense." 

The  bottom  velocity  amplitude  A  will  be  equated  to  U  in  the  analysis  of 
the  field  test  results.  The  bottom  velocity  amplitude,  U  ,  and  the  peak 
frequency, co  ,  are  tabulated  below  for  the  period  between  2400  hours 
18  February  and  0300  hours  21  February.  This  period  corresponds  to  the 
duration  of  a  mild  storm  and  constitutes  the  total  active  scouring  time 
between  13  February  and  28  February. 
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Date 

Time 

Bottom  Velocity  Amplitude 
um(in/sec) 

Peak  Frequency 

00  (cps) 

Depth  of*  Water 

4o  ft  60  ft 

d- =  0-221  mm  a  =  0.360  mm 

Depth  of  Water 

40  ft  60  ft 

18  Feb. 

2400 

3-934 

2.161 

0.1750 

0 . 1667 

19  Feb. 

0300 

5.169 

3.165 

0.1538 

0.1500 

0600 

5.566 

3.878 

0.1250 

0.1250 

0900 

9-584 

6.720 

0.1250 

0.1250 

1200 

9-6o4 

7.2o4 

O.II67 

0.1167 

1530 

8.304 

5.934 

O.U67 

0.1083 

1800 

6.744 

5.098 

O.U67 

0.1167 

2100 

5-546 

4.206 

O.I167 

0.1167 

2400 

5.565 

4.231 

O.H67 

0.1167 

20  Feb. 

0300 

5-011 

3.820 

O.II67 

0.1167 

0600 

4.583 

3-420 

0.1250 

0.1167 

0900 

4.002 

3-004 

O.II67 

0.1167 

1200 

3-241 

2.424 

O.II67 

0.1167 

1500 

3-011 

2-237 

O.I167 

0.1167 

1800 

2-437 

1-741 

0.1250 

0.1250 

2100 

2.147 

1.516 

0.1250 

0.1250 

2400 

1-957 

1.293 

0.1333 

0.1250 

21  Feb. 

0300 

i.943 

1.220 

0.1417 

0.1333 

Recorded  Settlement  -  The  position  of  each  mine  on  both  13  February 
and  28  February  are  as  listed  below: 
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for  t>ed-load  transport.  The  RHS  of  equations  ( 19a )  and  ( 19b )  is  equal 
to  the  RHS  of  equation  (10a)  or  (lla)  which  is  graphically  presented  as 
the  abscissa  on  Figure  13-  Hence,  upon  evaluating  the  LHS  of  equations 

ys 

(19a)  or  (l9b)  the  settlement,  -g—  ,  can  be  determined  from  Figure  13 
for  a  given  angle  of  repose,  4s  . 
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The  variation  in  the  bottom  velocity,  U  ,  during  the  period  of  the 
storm  is  shown  graphically  in  Figure  16  for  the  40-ft  water  depth  and 
in  Figure  17  for  the  60-ft  water  depth.  Linear  variation  of  U  with 
time,  t,  is  a  reasonable  representation  of  the  experimental  points  as 
shown  in  Figures  l6  and  17.  The  bottom  velocity,  U  ,  and  F  will  be 
represented  mathematically  by  these  straight  lines  in  order  that  the 
LHS  of  equations  (19a)  and  (l9b)  can  be  evaluated.  Since  the  sediment 
Froude  number,  F  =  U  Jyj  (8-1)  gd  ,  the  integrand  of  the  LHS  of 
equations  (19a)  and  (l9b)  can  be  written  as  F  and  F  ,  respectively. 

For  each  interval  on  Figures  1 6  and  17,  the  sediment  Froude  number 

F  =  a  +  bt 

in  which  a  is  a  dimensionless  constant  and  b  is  a  constant  having 
the  dimensions  of  sec  1.  The  values  of  a  and  b  for  their  various 
time  intervals  are  tabulated  below. 


From 

To 

a 

b 

(sec 

Time 

Interval, 

(sec) 

•ft  Depth: 

2400  hrs . 

18  Feb. 

0600  hrs.  19  Feb. 

1.69 

3.29(l0-'i) 

2.16(10^) 

0600 

19  Feb. 

0900 

19  Feb. 

2.40 

i.6o(io_lt) 

1.08(10^) 

0900 

19  Feb. 

1200 

19  Feb. 

4.l4 

0 

1.08(10^) 

1200 

19  Feb. 

2118 

19  Feb. 

4.i4 

-5.20(10”'’) 

3.35(104) 

2118 

19  Feb. 

2400 

19  Feb. 

2.40 

0 

9.72(io3) 

2400 

19  Feb. 

2100 

20  Feb. 

2.40 

-2.08(10"5) 

7.56(10^) 
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From 

To 

a 

b 

(sec 

Time 

Interval, 

(sec) 

-ft  Depth: 

2400  hrs. 

18  Feb. 

0600  hrs. 

19  Feb. 

0.75 

2.58(10-5) 

2.l6(l04) 

0600 

19  Feb. 

0900 

19  Feb. 

1.31 

8.84(10'5 ) 

1.08(104) 

0900 

19  Feb. 

1200 

19  Feb. 

2.27 

1.50(H)'5) 

1.08(104) 

1200 

19  Feb. 

20 40 

19  Feb. 

2.44 

-3.25(10'5) 

3.12(104) 

20 40 

19  Feb. 

2400 

19  Feb. 

1.42 

0 

1.20(l04) 

2400 

19  Feb. 

2100 

20  Feb. 

1.42 

-1.26(H)'5) 

7.56(l04) 

The  LHS  of  equations  (19&)  and  (19b)  can  he  readily  integrated  once 
the  decision  is  made  as  to  the  mode  of  transport.  The  criterion  for  the 
mode  of  transport  is  equation  (13),  for  which 


F  =  5-04  (  |  ) 

d  1/14 

at  transition.  For  values  of  F  below  ^.0h  (  —  )  ,  the  mode  of 

d  1/14 

transport  is  bed  load;  whereas  for  values  of  F  greater  than  5*04  (— ) 
the  mode  of  transport  is  suspended  load.  For  mines  A  and  B ,  which  were 
placed  in  the  60-ft  depth  the  mode  of  transport  throughout  the  storm  was 
bed  load  as  F  -  2.44.  For  mine  A,  which  was  placed  in  the  40-ft  depth, 
the  mode  of  transport  was  bed  load  up  until  0654  hours  on  19  February, 
suspended  load  from  0654  hours  to  1836  hours  19  February,  and  bed  load 
from  1836  hours  19  February  throughout  the  remainder  of  the  storm.  The 
settlement  of  mines  A,  B,  and  C  can  now  be  predicted  upon  integration  of 


the  appropriate  equations. 
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Mine  A  had  an  initial  burial  on  13  February  of  -g—  =  0 . 46o . 

Upon  considering  the  range  of  applicability  of  the  settlement  functions 

(Figure  15),  mine  A  falls  into  the  "no  settlement"  zone  for  the  highest 

values  of  F.  Since  mine  A  would  have  been  part  of  a  ripple  system  during 

the  period  of  highest  velocities  and,  since  the  velocities  during  the 

remainder  of  the  storm  were  quite  low,  no  scour  would  be  expected.  On 

y 

18  February  the  settlement  of  mine  A  was  -g—  =  0.460,  or  no  change. 

Mine  A  was  evidently  part  of  a  ripple  system  during  the  time  between 
13  February  and  28  February. 

Mine  B  ,  however,  did  experience  scour  during  the  storm.  It  was 

ysl 

repositioned  to  its  prone  position  (  -g—  =  0  )  on  13  February  and  on 

5S 

28  February  had  settled  to  -g—  =  0.405.  Upon  referring  to  Figure  15 
ysl 

for  d  =  O.360  mm,  rg^  =  0,  and  the  maximum  value  of  F  =  2.44,  mine  B 
is  seen  to  be  in  the  settlement  region  for  bed-load  transport.  The 
cumulative  value  of  equation  (l9h)  is  obtained  by  integrating  the  LHS 
between  the  limits  of  applicability  of  each  straight-line  segment  for  F 
and  then  summing  the  values  of  the  integrals  from  the  beginning  of  the 
storm.  The  cumulative  value  of  equation  (l9b)  and  the  corresponding 

ys 

values  of  -g—  are  tabulated  below.  The  predicted  and  observed  settle¬ 
ment  for  mines  A  and  B  are  also  shown  in  Figure  17- 

Mine  C  ,  which  was  placed  in  40-ft  of  water,  also  underwent  scour 
ys 

as  -g—  was  0.205  on  13  February  and  0.955  on  28  February.  Mine  C 

ysl 

lies  in  the  settlement  region  of  Figure  15  as  -g—  =  0.205,  and  d  =  0.221  nm. 
The  maximum  value  of  F  is  4.l4,  meaning  that  the  mode  of  transport  was 
suspended  load  during  the  most  severe  wave  action  of  the  storm.  Suspended- 
load  transport  prevails  for  f>2.88,  corresponding  to  the  period  between 
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0654  hours  and  1836  hours  19  February.  This  period  is  indicated  on 
Figure  l6.  Inasmuch  as  mine  C  was  initially  buried  on  13  February, 
and  the  mode  of  transport  initially  was  bed  load,  the  constant  of 

O.63 

ysl 

integration,  5-2  (  ~g—  )  >  must  be  included  in  the  LHS  of  equation  (19b). 

The  value  of  the  constant  of  integration  is  I.89.  Between  2400  hours 
18  February  and  0654  hours  19  February,  equation  (19b)  is  utilized.  From 
0654  hours  to  1836  hours  19  February  the  suspended-load  equation  (19a)  is 
integrated  and  the  values  of  the  integrals  cumulatively  added  to  the 
value  of  equation  (l9b),  obtained  previously.  From  1836  hours  19  February 
throughout  the  remainder  of  the  storm  equation  (19b)  is  utilized  as  bed¬ 


load  transport  prevails .  The  cumulative  values  of  the  LHS  of  equation 

ys 

(19a)  or  (19b)  and  the  associated  values  of  -g—  are  tabulated  below. 
The  predicted  and  observed  settlement  of  mine  C  are  shown  in  Figure  l6. 


Time 


Cumulative 
Value  of  LHS 
of  Equation 
(19a)  or  (19b) 


Predicted 

Settlement 

^s_ 

D 


Recorded 

Settlement 

^s_ 

D 


Mine  A: 


1200  hrs . 
1200 

Mine  B: 

13  Feb. 

28  Feb. 

Mine  A  is  pre¬ 
dicted  to  have 
been  part  of  a 
ripple  system 

O  O 
4=-  V 
&& 

1200  hrs. 

13  Feb. 

0 

0 

0 

2400 

18  Feb. 

0 

0 

- 

0600 

19  Feb. 

0.001 

0.040 

- 

0900 

19  Feb. 

0.079 

0.120 

- 

1200 

19  Feb. 

0.587 

0.240 

- 

20 40 

19  Feb. 

1.032 

0.325 

- 

2400 

19  Feb. 

1.033 

0.325 

- 

2100 

20  Feb. 

1.039 

0.330 

- 

1200 

28  Feb. 

1.039 

0.330 

0.405 
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Time 


Cumulative 
Value  of  IHS 
of  Equation 


(19a)  or  (19b) 


Predicted 

Settlement 

D 


Recorded 

Settlement 

^s_ 

D 


Mine  C: 


1200  hrs. 

13  Feb. 

- 

0.205 (given) 

0.205 

2400 

18  Feb. 

- 

0.205 (given) 

- 

0600 

19  Feb. 

2.13 

0.435 

- 

0654 

19  Feb. 

2.42 

0.470 

- 

0900 

19  Feb. 

4.36 

0.595 

- 

1200 

19  Feb. 

8.l4 

0.755 

- 

1836 

19  Feb. 

14.24 

0.935 

- 

2118 

19  Feb. 

15-04 

0.950 

- 

2400 

19  Feb. 

15.36 

0.955 

- 

2100 

20  Feb. 

15.74 

O.965 

- 

1200 

28  Feb. 

15-74 

0.965 

O.955 

The  agreement  between  predicted  and  observed  settlement  for  mines  A, 
B,  and  C  is  quite  satisfactory.  The  settlement  functions  derived  from 
the  model  tests  appear  to  be  suitable  for  prediction  of  prototype  mine 
settlement. 
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Predicted  and  Observed  Settlement  of  a  Mark-39  Mine  in 


Predicted  Seitleue; 


Conpariscn  of  Predicted  and  Observed  Settlement  of  Two  Mark-36 
Mines  in  60-ft  of  Water. 


CONFIDENTIAL 


CONCLUSIONS 

The  following  conclusions  were  reached  from  the  model  study  and  the 
associated  analyses. 

(1)  Mine  burial  occurs  as  the  result  of  (a)  scour  around  the  mine, 

(b)  settlement  of  the  mine  into  the  scour  hole,  and  (c)  subsequent  refilling 
of  the  scour  hole  between  storms. 

(2)  Scour  around  the  mine  occurs  either  by  the  material  being  removed 
as  bed  load  or  as  suspended  load.  The  settlement  functions  have  been  formu¬ 
lated  as  equation  (10a)  for  removal  as  suspended  load  and  as  equation  (lla) 
for  removal  as  bed  load.  The  limit  between  these  functions  is  given  by 
equation  (13 ) • 

(3)  As  settlement  occurs,  turning  also  occurs  such  that  the  mine  axis 
coincides  with  the  direction  of  the  wave  crest.  The  initial  orientation  has 
a  minor  effect  upon  settlement  for  initial  angles  of  60  degrees  or  less. 

(A)  Initial  burial  was  found  to  have  no  effect  upon  settlement  if  the 
material  is  removed  in  suspension.  Initial  burial  was  found  to  have  an 
effect  upon  settlement  if  the  material  is  removed  as  bed  load.  The  initial 
burial  constitutes  an  initial  condition  which  has  been  incorporated  into  the 
settlement  function,  equation  (lla) . 

(5)  Under  certain  flow  conditions  the  mine  becomes  part  of  the  ripple 
system  on  the  bed.  No  settlement  occurs  when  the  mine  is  a  pseudo  ripple. 
The  region  of  "no  settlement"  is  poorly  defined.  An  attempt  to  bound  this 
region  is  shown  in  Figure  15  for  an  l8-in.  diameter  mine.  Additional  ex¬ 
perimental  studies  will  be  required  in  order  to  more  precisely  establish  the 
bounds. 
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APPENDIX  A 


1 


SCOUR 


M.  R,  Carstens  and  C.  S.  Martin 


Scour  is  the  excavation  and  removal  of  bed  material  by  water  in 
motion.  A  general  mathematical  expression  can  be  formulated  from  the 
principle  of  conservation  of  mass  which  encompasses  all  scouring  situations 


Q  -  Q  . 

so  si 


d¥ 

dt 


(1) 


In  equation  (l),  Qsq  is  the  rate  that  sediment  is  being  removed  from 

the  scour  hole,  is  the  rate  that  sediment  is  being  transported  into 

the  scour  hole,  V  is  volume  of  the  scour  hole,  and  dv/dt  is  the  rate  of 

change  of  V.  A  uniform  stream  such  as  a  canal  section  is  stable  when 

the  LH3  of  equation  (l)  is  zero.  A  canal  is  scouring  when  the  capacity 

of  the  stream  to  excavate  and  remove  sediment  is  greater  than  the  rate 

of  sediment  inflow.  In  such  a  case  Q  >  Q  .  for  which  dv/dt  >  0.  The 

so  si 

increase  in  volume  -V- of  a  reach  is  accomplished  by  a  degradation  of  the 
bed  and/or  an  increase  in  the  canal  width. 

Localized  scour  occurs  in  the  vicinity  of  obstructions  placed  in 
the  flow.  The  increase  in  velocity  adjacent  to  an  obstruction  is 
accompanied  by  an  increased  capacity  to  carry  sediment  as  compared  with 
unobstructed  areas  of  the  bed.  In  many  situations  a  scour  hole  adjacent 
to  the  obstruction  will  occur  as  a  result  of  this  localized  increase  of 


2 


capacity  to  excavate  and  remove  bed  material. 

The  analysis  of  localized  scour  can  be  simplified  by  considering 
only  the  case  in  which  sediment  transport  into  the  scour  hole  is 
negligible.  In  this  case  the  sediment-transport  rate  out,  Qso>  is  equal 
to  the  rate  of  change  of  scour-hole  volume.  It  follows  from  equation  (l) 
that 

*  -  C  Qs0  dt  (2) 

J  o 

The  analysis  of  scour  can  be  systematized  by  considering  separately 
(a)  the  beginning,  (b)  the  active  phase,  and  (c)  the  termination.  Both 
the  beginning  and  termination  of  scour  are  dependent  upon  the  condition 
of  incipient  motion  of  the  bed  particles.  The  active  phase  of  scour  is 
concerned  with  the  scour-hole  geometry  and  the  sediment-transport  function. 

Scour-Hole  Geometry 

The  mathematical  description  of  scour-hole  geometry  in  uncon¬ 
solidated  materials  is  greatly  simplified  because  the  angle  between  the 
sides  of  the  hole  and  the  horizontal  is  the  angle  of  repose,  O,  of  the 

bed  material.  Hence  the  scour-hole  volume  is  proportional  to  the  depth 

2  3 

squared,  y  ,  in  a  two-dimensional  case  and  to  the  depth  cubed,  yg,  in  a 

three-dimensional  case.  From  equation,  (l)  with  no  sediment  inflow 


q  <*  y 
Mso  s 


% 

dt 


(3) 
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for  a  two-dimensional  hole  and 


o  o<  y2  dVs 

so  s  d r 


(4) 


for  a  three-dimensional  hole. 

An  excellent  example  which  illustrates  the  comparative  ease  of 
formulating  a  scour-hole  geometry  function  has  been  reported  by  Rouse 
(l).  In  Rouse's  experiment  a  submerged  two-dimensional  jet  was 
directed  down  a  vertical  wall  into  an  originally  level  sand  bed.  As 
would  be  expected  the  measured  scour  holes  were  geometrically  similar 
at  successive  times.  The  scour  hole  at  any  time  had  essentially  the 
configuration  shown  in  Figure  1. 


From  Figure  1  and  equation  (3) 

d¥  _  (l  +  1 V  1  \  dys 

%o  dt  \  {2/\tan  /  ys  dt 


(5) 
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in  which  the  scour  hole  extends  to  the  crest  of  the  dune. 

A  second  example  of  two-dimensional  scour  is  reported  by  Laursen 
(2).  Laursen1 s  experiments  were  very  similar  to  Rouse's  except 
that  the  two-dimensional  jet  was  directed  over  the  top  of  an  originally 
level  bed.  The  scour  hole  at  any  time  had  essentially  the  configuration 
shown  in  Figure  2. 


in  which  the  scour  hole  extends  to  the  crest  of  the  dune. 

An  example  of  three-dimensional  scour  is  reported  by  Ahmad  (3). 

In  one  experiment,  Ahmad  placed  a  vertical  plane  wall  in  otherwise  uniform 
open-channel  flow.  The  scour  hole  in  the  movable  bed  can  be  closely 
approximated  by  an  inverted  cone  with  the  apex  of  the  cone  at  the  end  of 
the  wall  as  in  Figure  3. 
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V 


a.  Top  View 


Figure  3.  Scour-Hole  of  Ahmad's  Experiment 
From  Figure  3  and  equation  (4) 

0  =  d£  =  — B _ 2  %  (7) 

vso  dt  tan2i>  's  dt 

Sediment-Transport  Functions 

General 

In  order  to  apply  equation  (2)  for  prediction  of  scour-hole 
development  with  time,  the  sediment- transport  rate  Qgo  must  be  formulated 
as  a  function  of  fluid  properties,  flow  characteristics,  sediment 
properties,  and  scour-hole  development.  The  sediment- transport  functions 
which  have  been  formulated  for  uniform  streams  are  not  suitable  as  trans¬ 
port  functions  for  localized  scour  for  two  reasons.  First,  the  uniform- 
stream  transport  functions  are  based  upon  the  average  shearing  stress 
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between  the  bed  and  the  fluid.  Inasmuch  as  localized  scour  occurs  in  non- 
uniform  flow  regions,  the  spatially-varying  shear  stress  at  the  surface 
of  the  bed  is  difficult  to  predict.  Second,  localized  scour  will  occur 
where  the  fluid  has  been  accelerated  and  the  boundary-layer  thickness  is 
negligible  in  contrast  to  uniform  stream  flows. 

With  a  negligible  boundary  layer,  the  transport  function  can  be 
formulated  qualitatively  without  reference  to  the  fluid  viscosity. 

The  drag  force,  F  ,  of  the  fluid  on  the  particle 

Fd  «  d2  pV2  (8) 

in  which  p  is  the  fluid  density,  and  V  is  the  fluid  velocity  at  the  top 
of  the  particle. 

The  resistance  of  particles  to  motion  can  readily  be  determined 
by  observation  of  the  angle  of  repose  J  of  a  submerged  pile  of  the  bed 
material  in  quiescent  water.  The  force  of  gravity  on  a  particle  down  the 
face  of  the  pile  is  resisted  by  the  particle-to-particle  reactions  of 
adjacent  particles.  Thus  for  particles  of  diameter,  d, 

c<  (ys~y)d^  sin  <J>  (9) 

— i 

in  which  F^  is  the  resisting  force,  y^  is  the  specific  weight  of  the. 

particles  and  y  is  the  specific  weight  of  the  water. 

The  local  sediment  transport  rate  AQg  will  be  a  function  of  the 

— ►  — ^ 

ratio  of  the  motivating  force,  F^,  to  the  resisting  force,  F^.  Since 
— 1 

the  motivating  force  is  proportional  to  the  fluid  inertial  force, 
equation  (8),  and  the  resisting  force  is  proportional  to  the  gravity 


force  on  the  particle,  equation  (9),  the  ratio  of  the  two  forces  is  a 
sediment  Froude  number 
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Thus 


(10) 


The  sediment  Froude  number,  F,  and  local  sediment  transport  rate,  AQs>  will 

vary  over  the  surface  of  the  scour  hole.  The  greatest  rate  of  transport 

will  occur  where  the  fluid  velocity  is  greatest  which  consequently  will  be 

the  position  of  greatest  scour  depth.  Since  the  transport  capacity  will 

decrease  away  from  the  position  of  greatest  depth,  yr,  much  of  the  material 

scoured  at  the  bottom  will  deposit  on  the  sides  of  the  scour  hole  and 

slide  down  toward  the  bottom.  Hence  the  net  transport  rate  Qsq  out  of 

the  hole  is  the  transport  out  of  the  periphery  of  the  hole.  As  the  scour 

hole  deepens,  the  lateral  limit  of  the  hole  moves  further  from  the  flow 

disturbance.  As  a  result  the  transport  capacity  at  the  edge  of  the  hole 

decreases  as  y  increases, 
s 

Thus 


Qs0  =  f  (  F  ,  yc,  geometry  ) 


(11) 


in  which  the  variable,  geometry,  pertains  to  the  geometry  causing  the  flow 
disturbance, 

A  different  process  of  removal  might  be  anticipated  either  when  F 
is  very  large  or  at  the  beginning  of  scour  when  the  hole  is  small.  As 
the  water  flows  over  the  bed  material,  some  material  may  be  transported 
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vertically  as  suspended  material.  If  this  material  leaves  the  scour  hole 
without  being  in  contact  with  the  bed  at  the  rim  of  the  hole,  the  functional 
relation,  equation  (11),  would  still  be  applicable  but  the  function,  per 
se .  would  be  different.  In  other  words,  different  sediment-transport 
functions  are  anticipated  for  removal  by  suspended  load  through  the  top 
than  for  removal  by  bed  load  over  the  rim  of  the  scour  hole. 

Analysis  of  Scour  Experiments 

Sediment-transport  functions  can  be  formulated  from  the  jet-scour 
experiments  of  Laursen  (2).  The  scour-hole  geometry  is  shown  schematically 
in  Figure  2.  The  results  are  presented  such  that  the  depth  of  scour,  y^, 
as  a  function  of  time,  t,  can  be  readily  determined  for  each  of  the  17 
runs.  The  jet  velocity  was  varied  from  run  to  run.  The  side-elevation 
view  of  the  scour  hole  was  obtained  photographically  at  successive  times 
during  a  run.  A  set  of  runs  was  repeated  for  each  of  three  sand  sizes. 

Using  Laursen' s  results  and  equation  (6),  the  sediment-transport  rate  was 
calculated  by  the  writers.  The  calculated  results  are  presented  in 
Figures  4  and  5. 

In  analyzing  Laursen1 s  results,  different  modes  of  transport  were 
found.  During  the  earlier  portions  of  15  runs  and  during  the  entire  run 
with  the  finest  sand  and  the  highest  velocity,  the  sediment-transport  rate 
was  found  to  be  independent  of  scour-hole  depth,  y  j  whereas  during  the 
later  portions  of  the  15  runs  and  during  the  entire  run  with  the  coarsest 
sand  and  lowest  velocity,  the  sediment-transport  rate  was  found  to  be 
dependent  upon  the  value  of  yg. 

The  writers'  interpretation  of  this  finding  is  as  follows.  In 
the  vicinity  where  the  jet  impinges  on  the  bed,  the  bed  material  is 
excavated  and  thrown  into  suspension.  This  suspended  material  is  trans- 
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Figure  5-  Sediment-Transport  Function  for  Laursen's  Bed-Load  Mode 
of  Transport. 
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ported  away  from  the  area  of  jet  impingement  while  settling  back  to  the 
bed.  During  the  earlier  portion  of  the  runs,  the  scour  hole  is  small 
enough  for  the  material  to  be  carried  out  of  the  hole  without  contacting 
the  bed.  In  this  situation  the  transport  rate,  q  ,  is  directly 
proportional  to  the  rate  at  which  the  jet  can  excavate  material  and  is 
independent  of  yg.  The  excavation  function  can  be  approximated  from 
Figure  4  as 

—  =  7  (10'4)  F  (12) 

q  o 

in  which  Fq  is  the  sediment  Froude  number  based  upon  'the  jet  velocity, 

V  . 
o 

As  the  scour  hole  deepens  and  the  dune  crest  moves  further  from  the 
jet-impingement  area,  the  suspended  material  settles  to  the  bed  before 
being  transported  over  the  crest.  In  this  case  the  transport  out  of  the 
hole  is  bed-load  transport  over  the  crest.  Since  the  jet  diffuses  in 
traveling  over  the  scour-hole  wall  the  bed-load  transport  capacity 
decreases  as  the  dune  crest  moves  further  from  the  jet-impingement  area. 

The  bed-load  transport  function  can  be  approximated  from  Figure  5  as 

C(  ¥  T  ( £  T  «> 

in  which  d  is  the  mean  diameter  of  the  bed  particles. 

A  transition  region  exists  between  the  suspended-load  transport 
and  bed-load  transport  conditions.  In  the  transition  region,  part  of 
the  excavated  material  is  transported  out  of  the  hole  as  suspended  load 
and  part  as  bed  load.  The  transition  region  is  surprisingly  narrow. 
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Ignoring  the  transition  region  completely,  the  limits  of  equation  (12)  are 

0  <  y  /b  <  0.78  (d/b)1/5  F  6/5 
s  o 

and  the  limit  of  equation  (13)  is 

ys/b  >0.78  (d/b)1//5  Fq 

Utilizing  equations  (12)  and  (13),  equation  (6)  can  be  integrated 
to  obtain  the  scour-hole  depth  as  a  function  of  time.  Assuming  yg  is  zero 
when  t  is  zero,  the  suspended  load  scour  function  is  of  the  form 

ys  =  k2tl/2  (14) 

If  yg  exceeds  the  limits  stated  above  then  both  the  suspended-load  trans¬ 
port  and  the  bed-load  transport  functions  must  be  employed  in  the 
integration  from  which  a  scour  function  of  the  form 

ys  =  [V  "  k3  ]  ^  {15) 

is  obtained.  The  point  to  be  emphasized  here  is  that  the  scour  depth  at 
any  time  is  dependent  upon  the  history.  In  Laursen's  experiments  the 
change  from  a  suspended-load  mode  of  removal  to  a  bed-load  mode  of 
removal  results  in  a  complicated  scour-time  function,  equation  (15). 

Sediment-transport  functions  can  also  be  formulated  from  the  jet 
scour  experiments  of  Rouse  (1).  The  scour-hole  geometry  is  shown 
schematically  in  Figure  1.  Rouse's  experiments  differ  from  Laursen's 
in  the  following  respects:  (a)  the  two-dimensional  jet  was  directed 
vertically  downward  onto  the  bed,  (b)  the  progress  of  scour  was  recorded 
by  wax-pencil  marking  on  the  glass  wall  of  the  flume,  and  (c)  one-half 
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scale  tests  were  performed  in  which  all  pertinent  geometric  dimensions 
were  halved  except  for  the  sand  size.  Rouse  observed  two  distinct  flow 
patterns  which  he  called  minimum  jet  deflection  and  maximum  jet  deflection. 
With  minimum  jet  deflection  the  main  stream  followed  up  the  scour-hole  slope. 
With  maximum  jet  deflection,  the  main  stream  separated  from  the  scour  hole 
forming  a  large  eddy  with  return  flow  down  the  wall  parallel  to  the  jet. 

Using  Rouse's  results  and  equation  (5)  the  sediment-transport  rate 
was  calculated  by  the  writers.  The  calculated  results  are  presented  in 
Figure  6  for  minimum  jet  deflection  and  in  Figure  7  for  maximum  jet 
deflection.  The  sediment- transport  function  for  minimum  jet  deflection 
can  be  approximated,  Figure  6,  as  follows 

rz2-5 (i°'6) fo8(?)  5  (16) 

Similarly  the  sediment- transport  function  for  maximum  jet  deflection  can 
be  approximated,  Figure  7,  as  follows 


—  =  (10-2)  F 

q 


o8(^) 


-10 


(17) 


The  condition  of  maximum  jet  deflection  occurred  initially  and  then 
switched  to  the  condition  of  minimum  jet  deflection  when 


^  =  0.064  FQ3(f)  (18) 

Both  the  initial  sediment-transport  function,  equation  (17),  and 
the  later  sediment- transport  function,  equation  (16),  appear  to  be  bed¬ 
load  functions.  This  conclusion  is  reached  by  comparison  with  the 


Figure  6.  Sediment-Transport  Function  for  Rouse's  Minimum  Jet 
Deflection. 
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sediment-tran sport  function  obtained  from  Laursen's  experiments.  In  other 
words,  equations  (l6)  and  (17)  closely  resemble  equation  (13)  but  are 
quite  dissimilar  to  the  suspended-load  function,  equation  (12).  The 
reason  is  that  during  the  early  stages  of  a  run  the  material  excavated 
by  the  vertically  plunging  jet  was  carried  into  the  recirculating  eddy 
adjacent  to  the  wall.  This  suspended  material  either  recirculated  or 
settled  back  into  the  scour  hole  rather  than  being  carried  out  of  the  hole 
as  suspended  load.  In  contrast  to  Rouse's  experiment,  Laursen  prevented 
the  separation  of  the  jet  stream  from  the  scour  hole  by  placing  a  lip  on 
the  upper  edge  of  the  slot  through  which  the  jet  issued. 

Utilizing  equations  (l6)  and  (17),  equation  (5)  can  be  integrated 
to  obtain  the  scour-hole  depth  as  a  function  of  time.  Assuming  y^  is 
zero  when  t  is  zero, 


y 


s 


tl/i2 


(19) 


if  y  is  less  than  the  limit  obtained  from  equation  ( 18 ) .  If  y  is 
s  *  s 

greater  than  this  limit, 

ys  =  [k5  1  ”  k6  ]  V7  (20) 


in  which  the  k's  are  functions  of  F  and  d/b. 

o 

Results  of  Ahmad's  experiments  (3)  can  be  analyzed  in  a  similar 
manner  to  the  analysis  performed  on  Laursen's  and  Rouse's  data.  Ahmad's 
experiments  were  conducted  with  bed  material  of  0.250  mm  diameter.  The 
vertical  wall  extended  to  one-half  the  channel  width  in  all  five  runs. 

The  variation  from  run  to  run  consisted  of  varying  the  channel  width  from 
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1.0  ft  to  3.0  ft  in  0.5  ft  increments.  The  sediment  ,Froude  number, 

V  /  J ( s-l) gd  ,  has  been  calculated  using  a  velocity  obtained  by  dividing 
the  discharge  by  the  cross-sectional  area  at  the  constricted  section. 
Details  of  the  experimental  procedure  are  not  given  in  reference  (3) . 

The  results  are  graphically  presented  from  which  depth  of  scour,  y  ,  as 
a  function  of  time  can  be  obtained.  Using  these  results  the  transport 
parameter,  G^/Vy^,  was  computed.  The  variation  of  the  transport  parameter 
is  shown  in  Figure  8.  The  sediment  transport  function  can  be  approximated 
as 


3.5  (10  U)  F8  ys  3 


(21) 


3 

The  coefficient  of  equation  (21)  has  the  dimensions  of  ft  .  The  only 
suitable  length  with  which  to  make  the  coefficient  dimensionless  is  the 
sediment  diameter,  d.  In  dimensionless  form,  equation  (21)  is 


0.063  F8  (  —  ^ 
V  d  y 


(22) 


Utilizing  equation  (22),  equation  (7)  can  be  integrated  to  obtain 
the  scour-hole  depth  as  a  function  of  time.  Assuming  y  =  0  when  t  -  0 

Ys  =  k?  t  1/4  (23) 


In  a  recent  paper  by  Laursen  (4)  the  history  of  scour  is  presented  for  a 
situation  similar  to  that  of  Ahmad's  experiments.  The  history  was 
determined  experimentally  at  Colorado  State  University.  The  chronological 
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development  of  the  scour  hole  is  such  that  is  proportional  to  the  fifth 
root  of  t  rather  than  the  fourth  root  as  indicated  in  equation  (23). 
Obviously  more  experimental  studies  are  needed  for  the  precise  definition 
of  the  sediment  transport  functions  of  scour. 

Conclusions 

Two  very  carefully  executed  experiments  involving  scour  by  jets 
were  analyzed  in  order  to  determine  the  nature  of  sediment- transport 
functions  for  scouring  situations  in  which  the  boundary-layer  thickness 
is  negligible.  The  conclusions  from  this  analysis  are  as  follows: 

1,  The  sediment  Froude  number  is  a  satisfactory  similarity  parameter 
for  model  studies  of  scour.  The  pertinent  properties  of  the  sediment  are 
correctly  incorporated  in  F  as  evidenced  by  the  minimum-deflection  trans¬ 
port  function  obtained  from  Rouse's  experiment,  that  is,  equation  (16). 

On  the  other  hand,  the  inclusion  of  F  did  not  completely  include  the 
sediment  properties  in  the  bed-load  function,  equation  ( 13) ,  obtained 
from  Laursen's  experiments.  Even  here  the  large  exponent  of  F  in 
comparison  to  the  low  exponent  on  the  term  containing  sediment  diameter, 
d,  is  reassuring.  The  sediment  Froude  number  is  also  a  logical  parameter 
for  suspended-load  movement  inasmuch  as  the  major  forces  on  a  particle 
in  suspension  in  a  turbulent  fluid  are  the  gravity  force  and  the  fluid- 
inertia  force.  Equation  (12)  and  Figure  4  support  this  observation. 

2.  The  bed-load  transport  functions,  equations  (13)  and  (l6), 
appear  to  be  typical  of  most  scouring  situations  in  that  the  transport 
is  by  bed  load  over  the  rim  of  the  scour  hole.  The  bed-load  transport 
function,  equation  (22),  for  a  three-dimensional  scour  situation  has  the 
same  functional  relationship  as  that  for  a  two-dimensional  scour  situation, 
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confirming  the  inclusion  of  the  sediment  Froude  number  to  a  large  exponent. 

3.  The  fact  that  scour-hole  dimensions  at  any  time  are  dependent 
upon  history  can  lead  to  difficulties  in  analysis  of  experimental  results. 
Since  scour-hole  dimensions  are  proportional  to  the  integral  of  transport 
rate  with  time,  equation  (2),  interpretation  of  the  integral  can  be 
difficult.  For  example,  in  Laursen's  experiments  the  removal  was  initially 
as  suspended  load,  equation  (12),  and  later  as  bed-load,  equation  ( 13) . 

The  writers  avoided  this  difficulty  by  analyzing  the  rate  of  removal  which 
is  dependent  upon  prevailing  flow,  fluid,  and  geometric  conditions. 


in  which  ZF,,  is  the  summation  of  the  forces  parallel  to  the  bed  and  ZF 
is  summation  of  the  forces  perpendicular  to  the  bed.  The  parallel  forces 
are  the  drag  force  on  the  particle  and  the  component  of  the  gravity  force 
down  the  plane.  The  perpendicular  force  is  the  component  of  gravity  force. 


Thus 
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d  pV/2-  (y  -  y)  d  sir>  a 

tan  O  c<  - - - (25) 

(ys  ~  Y)a  cos  a 

in  which  a  positive  value  of  a  is  associated  with  flow  up  the  slope. 

Upon  rewriting  this  relationship  a  parameter  similar  to  the  sediment 
Froude  number  appears 

F  c<  V  =  C  (26) 

</(s-l)gd  [tan  $  cos  a  +  sin  a] 

This  parameter  should  have  a  constant  value,  C,  for  incipient  scour  for 
which  there  is  a  negligible  boundary  layer.  The  critical  scour  parameter, 
C,  will  be  computed  from  results  of  several  experiments  in  order  to  de¬ 
lineate  its  value  for  incipient  scour. 

White  (5)  determined  the  conditions  for  incipient  motion  of  bed 
materials  in  a  converging  flow  of  both  air  and  water.  With  a  converging 
flow  and  low  viscosity  fluids,  the  boundary  layer  is  quite  thin,  A 
direct  measure  of  the  relative  size  of  a  sand  grain  to  the  thickness  of 
a  turbulent  boundary  layer  is  the  shearing-velocity  Reynolds  number 

R  d. 

*  v  v 

in  which  is  the  critical  boundary  shear  stress,  p  is  the  density  of 
the  fluid,  and  V#  is  the  defined  shearing  velocity.  As  increases  the 
boundary-layer  thickness  becomes  smaller  in  comparison  with  the  sand  grain 
diameter,  d.  With  the  comparatively  thin  boundary  layer  and  hence  no 
measurable  velocity  variation  over  the  bed  particle,  the  critical  scour 
parameter,  C,  should  attain  a  constant  value.  White's  experimental  results 
and  computed  values  of  the  critical  scour  parameter,  C,  are  listed  in  Table 
1.  The  critical  scour  parameter,  C,  is  essentially  constant  for  R#  of  80 


and  greater.  The  highest  values  of  correspond  to  flows  of  a  localized- 
scour  nature.  The  minimum  value  of  C  of  2.5  appears  to  be  suitable  for 
design  against  incipient  scour. 

Ippen  and  Verma  (6)  determined  the  conditions  for  incipient  motion 
of  isolated  spherical  particles  resting  on  a  horizontal  sand  bed.  A 
steady  stream  of  water  passed  over  the  particles.  Although  Ippen  and 
Verma  do  not  list  any  values  of  the  natural  angle  of  repose,  |3>,  the 
critical  scour  parameter,  C,  can  still  be  computed  upon  using  results 
given  by  Eagleson  and  Dean  (7).  Eagleson  and  Dean  determined  the  angle 
of  repose  for  a  similar  bed  configuration.  Experimental  results  from 
Ippen  and  Verma  and  computed  values  of  C  are  listed  in  Table  2.  The 
critical  scour  parameter  essentially  decreases  as  the  shearing-velocity 
Reynolds  number  increases,  as  with  the  results  of  White.  The  lowest  value 
of  C  of  1.5  is,  of  course,  less  than  the  corresponding  lowest  value  of 
2,5  for  White.  This  difference  is  due  mainly  to  the  fact  that  the  pro¬ 
jected  flow  area  of  the  isolated  spheres  of  Ippen  and  Verma  is  consider¬ 
ably  larger  than  that  of  the  imbedded  sand  particles  of  White.  Nevertheless, 
the  value  of  C  of  1.5  for  the  isolated  spherical  particles  should  represent 
an  absolute  minimum  for  localized  scour. 

Incipient  motion  tests  for  isolated  spherical  particles  of  high 
density  were  conducted  at  Georgia  Tech.  The  spherical  particles  were 
steel  and  the  fluid  was  water.  The  experiment  was  conducted  in  a  glass- 
walled,  enclosed,  rectangular-shaped  conduit  which  has  l/l6-in  diameter 
piezometer  holes  drilled  through  a  smooth  floor.  The  steel  particles 
(5/l6-in  diameter)  were  placed  such  that  they  rested  on  the  piezometer 
holes.  From  the  known  geometry  the  natural  angle  of  repose  can  be  com¬ 
puted.  Incipient  motion  was  determined  by  setting  the  flume  at  an  angle 
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a  with  the  horizontal  and  then  increasing  the  velocity  in  the  enclosed 
section  until  the  ball  jumped  out  of  the  hole.  Experimental  values  for 
incipient  motion  and  computed  values  cf  C  are  presented  in  Table  2.  The 
results  here  agree  quite  well  with  those  of  Ippen  and  Verma. 

To  further  illustrate  the  significance  of  the  critical  scour 
parameter,  experimental  results  for  incipient  motion  of  bed  particles  in 
oscillatory  flow  will  be  analyzed.  Manohar  (8)  determined  conditions  for 
incipient  motion  by  oscillating  a  bed  of  particles  simple  harmonically 
under  an  otherwise  still  tank  of  water.  The  boundary  layer  was 
turbulent  for  the  results  to  be  analyzed  here.  Manohar  conducted  many 
tests  for  each  bed  material  by  varying  the  amplitude  and  the  frequency 
of  oscillation  of  the  bed.  He  found  that,  for  each  bed  material, 
incipient  motion  occurred  at  a  constant  value  of  the  maximum  velocity  of 
the  bed,  Um<  Since  Manohar  did  not  experimentally  determine  the  angle 
of  repose,  3>,  of  his  bed  materials,  the  critical  scour  parameter  cannot 
be  easily  determined.  However,  the  writers  feel  that,  after  surveying 
experimental  results  presented  in  the  literature,  the  angle  of  repose  of 
Manohar1 s  bed  materials  is  most  probably  between  30°  and  40°.  Values  of 
the  critical  scour  parameter  based  on  the  maximum  velocity  are  computed 
with  the  angle  of  repose  assumed  to  be  30°  and  40°.  The  critical  scour 
parameter  and  Manohar1 s  experimental  results  for  incipient  motion  are 
presented  in  Table  3.  The  maximum  velocity,  U  ,  for  each  sand  is  the 
average  of  many  runs. 

As  with  the  results  of  White,  the  critical  scour  parameter  de¬ 
creases  as  the  bed  particle  diameter  increases  indicating  that  the  larger 
particles  are  not  influenced  by  the  boundary  layer  nearly  as  much  as  are 
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the  smaller  ones.  For  the  two  largest  particles  the  critical  scour  parameter 
is  essentially  constant  and  should  represent  the  condition  for  incipient 
scour.  An  average  value  of  2,8  for  C  would  seem  suitable  for  design.  This 
value  is  only  slightly  higher  than  the  2.5  resolved  from  White's  data. 

The  value  of  C  would  be  somewhat  higher  from  Manohar's  data  because  the 
maximum  velocity  Um  was  used  in  calculating  C.  Manohar  defined  incipient 
motion  as  being  when  a  few  particles  moved  a  short  distance.  Obviously 
the  velocity  would  have  to  be  greater  than  a  threshold  velocity  in  order 
for  the  particles  to  move  a  short  distance.  If  this  threshold  velocity  is 
0.9  the  C  values  obtained  from  Manohar's  and  White's  experiment  coincide. 
This  analysis  for  incipient  motion  in  oscillatory  flow  again  demonstrates 
the  significance  of  the  critical  scour  parameter  for  flows  with  negligible 
boundary  layers. 

Beginning  and  Termination  of  Scour 

Scour  occurs  whenever  and  wherever  the  critical  scour  parameter  is 
exceeded.  Once  exceeded  a  scour  hole  is  formed.  Scour  terminates  when  the 
bed  particles  can  no  longer  be  carried  over  the  crest  of  the  dune  by  the 
fluid.  Thus  the  critical  scour  parameter  is  the  governing  criterion  for 
both  the  beginning  and  the  end  of  scour.  The  critical  scour  parameter 
for  beginning  of  scour  is  computed  with  a  =  0  since  the  bed  is  not  yet 
deformed  by  scour.  The  critical  scour  parameter  for  terminal  scour  is 
computed  with  a  =  ®  since  the  scour  hole  slope  is  that  of  the  natural  angle 
of  repose  of  the  bed  material.  In  either  case  the  local  velocity  rather 
than  a  mean  must  be  employed  in  the  evaluation  of  the  critical  scour 
parameter,  C,  by  means  of  equation  (26). 

Laursen  performed  some  additional  experiments  from  which  the  value 
of  C  can  be  approximately  determined  when  scour  ceases.  Laursen  determined 
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the  minimum  jet  velocity,  V  ,  which  would  transport  particles  up  the  slope 
of  the  scour  hole  depicted  in  Figure  2,  but  not  over  the  crest.  From  one 
of  his  graphs  the  critical  scour  parameter  can  be  resolved.  The  velocity 
at  the  crest,  V,  which  is  considerably  less  than  V  ,  can  be  approximated 
from  a  relationship  given  by  Albertson,  ai  aJ.*(9)  for  the  diffusion  of  a 
two-dimensional  jet  from  a  slot.  In  terms  of  Laursen's  variables 


2.28 


cos  3> 


(27) 


in  which  XL  is  the  horizontal  distance  from  the  slot  to  the  crest  of  the 
dune.  The  distance  X^/cos  <3>  is  the  total  assumed  distance  traveled  by 
the  water  along  the  bed  of  the  scour  hole.  This  distance  corresponds  to 
the  centerline  distance  from  slot  to  the  point  in  question  for  the  study 
of  Albertson,  et  al.  With  the  angle  of  repose  of  32.6°,  with  b  of  0.025 
ft,  and  with  the  specific  gravity  of  2.65,  the  critical  scour  parameter, 

C,  was  computed  by  using  equation  (27)  to  obtain  V  and  by  using  a  =  4>  in 
equation  (26).  The  computed  values  are  presented  in  Table  4-  The 
critical  scour  parameter  is  essentially  a  constant  from  Laursen's  tests 
for  the  largest  values  of  X^,  but  considerably  less  than  that  resolved 
from  White's  results.  The  lower  value  of  C  for  Laursen's  tests  is 
believed  to  be  due  to  the  difference  in  flow  patterns  of  Albertson's 
and  Laursen's  tests.  In  other  words,  equation  (27)  is  only  qualitatively 
correct  for  application  to  Laursen's  scour-hole  geometry.  In  any  event 
the  critical  scour  parameter  is  a  constant  for  the  termination  of  scour. 

An  exceptionally  neat  experiment  which  demonstrates  the  usefulness 
of  the  critical  scour  parameter,  C,  is  reported  by  Chabert  and  Engeldinger 
(10)  .  In  these  experiments  water  flowed  around  a  single  vertical  cylinder 
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and  over  a  movable  bed  of  3-mm  diameter  gravel.  Three  different  cylinders 
were  used  with  diameters  of  5,  10,  and  15  cm.  The  water  depth  in  the 
approach  channel  was  systematically  varied  being  either  10  cm,  20  cm,  or 
35  cm.  The  remaining  variable  was  velocity,  V^,  approaching  the  cylinder. 
For  each  depth  of  flow  a  number  of  runs  were  made  with  the  velocity  at  a 
selected  value,  ranging  from  0.3  to  1.2  m  per  sec.  The  terminal  depth 
of  scour,  y  ,  was  observed  for  each  of  the  75  runs.  The  results  are 
plotted  in  dimensionless  form  in  Figure  9. 

The  beginning  of  scour  is  shown  at  a  sediment  Froude  number  of 
1.25  in  Figure  9.  For  3-mm  gravel  a  reasonable  value  of  O  is  45  degrees. 

At  the  beginning  of  scour,  the  bed  is  undeformed,  that  is  a  is  zero.  In 
these  circumstances,  the  critical  scour  parameter,  equation  (26),  is  simply 

C  =  — — -  (28) 

J  (s-l)gd 

With  3-mm  gravel  and  in  an  accelerated  flow  zone,  the  boundary  layer  would 
be  insignificant.  In  this  situation  a  value  of  C  of  2.5  is  to  be  expected. 
Scour  will  begin  at  the  point  of  highest  velocity  which  would  occur  at 
about  the  side  tangent  points  on  the  cylinder.  The  velocity  at  these  points 
can  be  taken  as 


V  =  2V  (29) 

from  potential  theory.  Upon  substitution  of  C  =  2.5  and  equation  (29)  into 
equation  (28),  the  value  of  J (s-i)gd  at  which  scour  can  be  expected 

to  begin  is  1.25.  The  experimental  results  shown  in  Figure  9  are  in 
excellent  agreement  with  the  prediction  of  the  beginning  of  scour. 


Figure  9-  Terminal  Scour  Depth  for  Flow  Around  a  Vertical  Cylinder 
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When  the  entire  bed  goes  into  motion,  sediment  is  carried  into  the 
scour  hole  at  the  same  rate  as  sediment  is  removed.  Thus  after  the  entire 
bed  is  in  motion  no  further  increase  is  expected  in  the  terminal  scour  depth, 
y  Again  the  experimental  results  are  in  excellent  agreement  with  the 
predicted  value  of  /  J ( s-1 ) gd  being  2.5.  The  value  of  y  ^/D  is  essentially 

constant  for  values  of  /  J ( s-l) gd  greater  than  2.5. 

When  the  value  of  /  */(s-l)gd  is  in  the  range  from  1.25  to  2.5 
localized  scour  occurs  around  the  cylinder  with  the  bed  being  motionless 
away  from  the  cylinder. 

Summary 

The  critical  scour  parameter  is  the  governing  parameter  for  incipient 
and  terminal  scour  since  boundary  layer  effects  are  negligible  where 
localized  scour  occurs.  The  angle  of  repose  is  a  variable  for  incipient 
motion  criteria.  For  the  prevention  of  scour  the  critical  scour  parameter 
should  not  exceed  2.5  in  unidirectional  flow. 

An  engineering  example  of  the  use  of  the  critical  scour  parameter 
C  is  the  determination  of  the  size  of  bed  material  to  be  placed  around  a 
cylindrical  pile  which  is  subject  to  wave-induced  flow.  The  waves  will 
be  assumed  to  be  simple  harmonic  in  shape.  The  still-water  depth  is  h, 
the  wave  amplitude  from  trough  to  crest  is  H  and  the  wave  length  is  L. 

From  Lamb  (ll)  the  maximum  velocity  of  the  water  at  the  bed, 

II  _  _ aU _  (go) 

o  2  c  cosh  2jth/L 

in  which  c  is  the  phase  velocity  of  the  wave.  If  the  flow  in  the 

accelerated  flow  region  adjacent  to  the  pile  is  assumed  to  be  irrotational , 

the  maximum  velocity  at  the  bed,  U  =  2U  .  If  the  value  of  the  critical 
'  mo 
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scour  parameter  of  2.8  from  Manohar's  results  is  used 

oH/c  cosh  2ith/L  -  2  8  (31) 

</ ( s-l)gd  tan  <P 

for  incipient  scour.  For  shallow-water  waves,  for  which  h/L  is  small, 
c  =sji/gh  and  cosh  2nh/L  1.  Thus 

■  ■  ~ H  -  =  2.8  (32) 

J ( s-l)hd  tan  0 

for  shallow-water  waves.  Equation  (32)  is  a  simple  formulation  for  the 
prevention  of  scour  at  a  cylindrical  pile.  This  example  shows  the  ease 
with  which  the  critical  scour  parameter  criterion  may  be  applied  to  a 
particular  situation  for  which  the  maximum  local  velocity  is  known.  This 
example  also  shows  the  large  diameter  of  the  particles  required  in  order 
to  prevent  scour  around  the  pile.  If  H  =  8  ft,  h  =  40  ft,  O  =  45°,  and 
s  =  2.6,  the  minimum  diameter  of  placed  material  would  be  3.9  cm. 


tables 
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Table  1.  Incipient-Motion  for  Unidirectional  Flow 
over  a  Homogeneous  Bed  (White's  Data) 


Bed 

Material 

d  (mm) 

Fluid 

ys/t 

a 

<£> 

V 

(ft/sec) 

V*d 

V 

C 

Sand 

0.90 

Water 

2.6 

0° 

45° 

1.25 

33 

3.2 

Steel  Shot 

0.71 

Water 

7.9 

0° 

34.4° 

2.03 

35 

3.4 

Sand 

0.90 

Air 

2100 

0° 

45° 

37.5 

80 

2.7 

Sand 

5.6 

Water 

2.6 

-24.2° 

45° 

1.84 

360 

2.7 

Sand 

5.6 

Water 

2.6 

0° 

45° 

2.46 

480 

2.5 

Sand 

5.6 

Water 

2.6 

26.2° 

45° 

2.98 

590 

2.6 

Sand 

5.6 

Air 

2100 

0° 

45° 

88.4 

1280 

2.5 
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Table  2.  Incipient-Motion  Results  for  Isolated  Spherical  Particles 


Investigator 


Ippen  and 
Verma 


Georgia 
Tech  Tests 


Material 

d(mm) 

ys/y 

a 

4> 

V 

V«d 

V 

C 

Plastic 

2.00 

1.28 

0° 

33.6° 

0.35 

18 

1.8 

Plastic 

2.00 

1.28 

0° 

33.6° 

0.36 

21 

1.8 

Plastic 

3.17 

1.28 

0° 

23.2° 

0.33 

27 

1.7 

Plastic 

3.17 

1.28 

0° 

23.2° 

0.35 

28 

1.7 

Plastic 

2.00 

1.28 

0° 

49.7° 

0.51 

30 

1.9 

Plastic 

2.00 

1.28 

0° 

49.7° 

0.52 

33 

2.0 

Plastic 

2.00 

1.28 

0° 

49.7° 

0.53 

34 

2.0 

Glass 

3.17 

2.38 

0° 

23.2° 

0.55 

42 

1.3 

Plastic 

3.17 

1.28 

0° 

42.7° 

0.46 

50 

1.6 

Plastic 

3.17 

1.28 

0° 

42.79 

0.50 

52 

1.7 

Glass 

3.17 

2.38 

0° 

23.2° 

0.69 

53 

1.6 

Plastic 

3.17 

1.28 

0° 

42.7° 

0.51 

54 

1.7 

Glass 

4.00 

2.38 

0° 

18.0° 

0.52 

55 

1.3 

Glass 

4.00 

2.38 

0° 

18.0° 

0.69 

67 

1.6 

Glass 

3.17 

2.38 

0° 

42.7° 

1.01 

80 

1.6 

Glass 

3.17 

2.38 

0° 

42.7° 

0.97 

80 

1.5 

Glass 

3.17 

2.38 

0° 

42.7° 

1.04 

83 

1.6 

Glass 

4.00 

2.38 

0° 

38.3° 

0.99 

100 

1.5 

Glass 

4.00 

2.38 

0° 

38.3° 

1.02 

103 

1.5 

Glass 

4.00 

2.38 

0° 

38.3° 

1.04 

106 

1.5 

Steel 

7.93 

7.63 

-6° 

11.3° 

0.97 

- 

1.8 

Steel 

7.93 

7.63 

-2° 

11.3° 

1.25 

- 

1.7 

Steel 

7.93 

7.63 

0° 

11.3° 

1.29 

- 

1.5 

Steel 

7.93 

7.63 

2° 

11.3° 

1.40 

- 

1.5 

33 


Table  3.  Incipient  Motion  Results  for  Oscillatory  Flow  Over  a 
Homogeneous  Bed  (Manohar's  Data) 


Material 

d(mm) 

Speci fic 
Gravity  (s) 

U  (ft/sec) 
m 

C 

0>  =  30° 

4>  =  40° 

Glass  Beads 

0.610 

2.54 

1.000 

4.2 

3.5 

Sand 

0.786 

2.63 

1.072 

3.8 

3.2 

Sand 

1.006 

2. 60 

1.129 

3.6 

3.0 

Sand 

1.829 

2.60 

1.355 

3.2 

2.7 

Sand 

1.981 

2.63 

1.380 

3.1 

2.6 

34 


Table  4*  Incipient-Motion  for  Jet  Flow  over  a  Homogeneous  Bed 

(Laursen's  Data) 


d  (mm) 

XL  (in) 

V  (ft/sec) 

C 

0.240 

7.4 

0.27 

1.3 

10.9 

0.31 

1.5 

14.7 

0.30 

1.4 

18.5 

0.32 

1.5 

20.3 

0.34 

1.6 

25.0 

0.32 

1.5 

32.9 

0.32 

1.5 

41.4 

0.32 

1.5 

0.690 

4.8 

0.39 

1.1 

6.0 

0.43 

1.2 

6.9 

0.44 

1.2 

8.5 

0.42 

1.2 

11.4 

0.53 

1.5 

13.1 

0.49 

1.4 

14.7 

0.53 

1.5 

16.5 

0.48 

1.3 

22.3 

0.50 

1.4 

25.0 

0.59 

1.6 

28.5 

0.61 

1.7 

31.8 

0.50 

1*4 

36.1 

0.62 

1.7 

39.5 

0.61 

1.7 

1.600 

6.0 

0.72 

1.3 

10.7 

0.75 

1.4 

14.1 

0.81 

1.4 

15.0 

0.84 

1.5 

18. 1 

0.89 

1.6 

21.8 

0.84 

1.5 

25.0 

0.77 

1.4 

26.8 

0.91 

1.7 

32.9 

0.78 

1.4 

37.0 

0.87 

1.6 

39.6 

0.90 

1.6 

35 
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APPENDIX  B 
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A  THEORETICAL  AND  EXPERIMENTAL  INVESTIGATION 
OF  FLOW  UNDER  A  PARTIALLY -IMBEDDED  CYLINDER 


C.  S.  Martin 


Flow  of  a  fluid  under  a  circular  cylinder  which  is  half-imbedded 
in  sand  is  amenable  to  theoretical  analysis  as  long  as  the  critical 
gradient  at  the  bed  is  not  attained.  The  theoretical  solution  for  a 
partially-imbedded  cylinder  is  untenable  in  that  singular  points  occur 
at  points  of  imposed  nonorthogonality  of  flow  lines  and  potential  lines. 
In  the  following,  an  experimental  study  shows  the  behavior  of  the  bed  in 
the  vicinity  of  the  singular  points.  The  effect  of  flow  out  of  or  into 
a  sand  bed  on  the  angle  at  which  the  bed  makes  with  the  horizontal  is 
discussed.  Experimental  results  are  given  on  vertical  piezometric -head 
gradients  required  to  cause  boiling  and,  finally,  piping. 


Theoretical  Solution 

In  Figure  1  is  shown  a  circular  cylinder  partially  imbedded  in 
sand.  Water  flows  through  the  sand  from  right  to  left  under  the  cylinder 
by  virtue  of  the  water  level  h ^  being  larger  than  h^.  The  problem  is 
to  determine  the  variation  of  the  piezometric  head,  h,  within  the  sand. 
Flow  through  homogeneous  porous  media  can  be  represented  by  Darcy' s  Law; 
namely, 


V 

x 


V 

y  = 


K 


ah 

ay 


2 


in  which  V  and  V  are  the  x-  and  y-components  of  macroscopic  velocity 
x  y 

in  the  bed,  and  k  is  the  coefficient  of  permeability. 


Figure  1.  Sketch  of  Imbedded  Cylinder 


Satisfying  the  equation  of  continuity  of  fluid  flow  through  the  sand 
results  in  Laplace's  equation 

aih  +  &  =  o  (i) 

2  2 

3x  3y 

or,  if  4>  =  -  Kh, 

,  3^  =  n  (2) 

2  2 

3x^  3y 

The  two  constant  piezometric  heads,  h^  and  h^,  make  the  sand-water 
interfaces  on  each  side  of  the  cylinder  potential  lines.  These  are 
boundary  conditions.  The  imbedded  portion  of  the  cylinder  is  regarded 


as  a  streamline.  The  solution  to  this  boundary-value  problem  is 
approached  by  the  technique  of  conformal  mapping.  The  physical  plane 
is  the  z-plane  in  which 
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z  =  x  +  i  y 

The  complex  function 

w  =  <P  +  i  t|> 


in  which  iji  is  the  stream  function,  represents  what  is  known  as  the  w- 
plane.  The  problem  is  to  map  the  w-plane  into  the  z-plane,  a  procedure 


which  requires  two  intermediate  planes,  the  z' -plane  and  the  t-plane. 


These  four  planes,  and  the  corresponding  points,  are  shown  in  Figure  2. 
In  Table  1  are  listed  the  complex  values  of  each  point  for  the  four 
planes. 

Table  1.  List  of  Complex  Values  of  the  Four  Planes 


Point 

z-plane 

z’ -plane 

t-plane 

w-plane 

A 

-  00 

0 

00 

*1  +  1  *A 

B 

-  Va(2r0  “  a) 

»  +  i(-  9/2  to  0) 

1 

+  i  '1'1 

C 

-  i  3 

-i8/2 

0 

®c  +  1  *1 

D 

~\Ja( 2r  -  a) 
v  0 

w  +  i(-  0/2  to  0) 

-1 

®2  +  i  +! 

1 

E 

CO 

0 

—  CO 

*2  +  1  +E 

The  transformation  from  the  z-plane  to  the  z'-plane,  as  listed  by 
Kober  ( l) ,  is 


5 


o 


or 

w  =  —  (  O  -  4>  )  cosh  +  0.  +111), 

ic  l  z  l  i 


upon  joining  the  w-plane  and  the  t-plane.  The  solution  connecting  the 
z-plane  to  the  w-plane  is 


The  vertical  piezometric-head  gradient,  8h/8y,  is  determined  by 
differentiating  the  real  part  of  equation  (4).  The  vertical  piezometric- 
head  gradient  on  the  bed 


It  is  of  interest  to  know  the  magnitude  of  8h/8y  at  point  B,  since  the 
sand  might  tend  to  become  "quick"  there.  Unfortunately,  however,  the 
value  of  8h/8y  is  mathematically  indeterminate  at  B  and  at  D. 

Except  for  the  case  in  which  the  cylinder  is  half-imbedded, 
potential  lines  AB  and  DE  are  not  orthogonal  with  streamline  BCD  at 
intersections  B  and  D.  Points  B  and  D  are  called  singular  points.  When¬ 
ever  orthogonality  is  not  satisfied,  the  piezometric-head  gradient  will 
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have  either  of  two  values,  as  discussed  by  Casagrande  (2).  He  states 
that  at  points  where  boundary  flow  lines  (streamlines)  intersect  potential 
lines  (lines  of  constant  piezometric  head)  at  a  predetermined  angle,  the 
piezometric-head  gradient  is  either  zero  or  theoretically  infinite.  If, 
as  shown  in  Figure  3,  the  angle  p  defines  the  angle  of  intersection  of 
boundary  flowline  and  potential  line,  criteria  from  Casagrande  may  be 
stated  as  follows 


If 

p<  90°  , 

the 

piezometric-head  gradient 
zero. 

is 

If 

P=  90°  , 

the 

piezometric-head  gradient  is 
finite,  not  necessarily  zero 

If 

p>  90°  ,  o 

the 

piezometric-head  gradient 
theoretically  infinite. 

is 

Figure  3.  Boundary  Angle  between  Streamline 
and  Potential  Line 

Casagrande  further  points  out  that  the  condition  of  infinite  gradient 
(or  infinite  velocity)  is  the  cause  of  erosion  in  dams  at  points  for 
which  p  >  90°.  Furthermore,  large  velocities  will  render  Darcy's  law 
invalid  since  velocity-head  changes  cannot  be  neglected. 
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Figure  4.  Sketch  of  Angle  of  Intersection  of 
Streamline  and  Potential  Line  for 
Cylinder 

Casagrande's  criteria  will  be  used  with  reference  to  Figure  4,  If 
a/rQ  >  1  ,  p  <  90°  ;  and  if  a/rQ  <  1  ,  p  >  90°  .  At  point  B  the 
piezometric-head  gradient  is  zero  for  a /tq  >  1  ,  finite  for  a/r^  =  1  , 
and  infinite  for  a/rQ  <  1  .  For  the  half-imbedded  case,  that  is, 
a/rQ  =  1  ,  the  vertical  piezometric-head  gradient  at  the  bed  is  given 
by 

Qh  |  =  h2~  hI  (6) 

3y  'y=0 

At  point  B, 

ah  |  =  h2  ~  hi 

ay  B  nr 

o 

If  the  mathematical  limit  is  indeed  a  physical  limit,  then  move¬ 
ment  of  the  bed  material  at  point  B  could  be  expected  whenever  the 
cylinder  is  less  than  one-half  buried  regardless  of  the  magnitude  of 

-  hf  .  From  the  physical  standpoint  this  conclusion  seems  preposterous. 
Only  by  experiment  can  this  dilemma  be  resolved. 
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Experimental  Setup 

A  photograph  of  the  experimental  apparatus  used  is  shown  in 
Figure  5.  The  walls  and  bottom  are  made  of  1/2-in  plexigla^  and  are 
glued  and  screwed  together.  The  flow  passage  is  4  in  wide.  The  black 
circular  disc  in  the  center  is  the  end  of  cylinder  which  is  4-in  in 
diameter.  The  lower  flow  boundary  is  a  plastic  strip  which  is  bent  in 
the  shape  of  a  12-in  circular  arc.  The  bed  material  is  placed  between 
the  cylinder  and  the  arc  up  to  any  desired  level  on  the  cylinder.  Water 
enters  through  an  orifice  in  the  right-hand  wall  and  flows  through  the 
sand  bed  under  the  cylinder.  A  constant  piezometric  head,  h^,  is  main¬ 
tained  on  the  left-hand  side  by  means  of  a  sharp-crested  weir.  For  a 
given  rate  of  flow  through  the  sand  the  piezometric  head,  reaches  a 
level  of  equilibrium,  provided  that  h^  -  h^  is  below  that  value  at  which 
piping  occurs. 


Permeability  Tests 

Tests  were  conducted  to  determine  the  coefficient  of  permeability, 
k,  of  each  bed  material.  The  physical  characteristics  of  the  two  bed 
materials  used  are  listed  below 


Material 


a  immJ 


Geometric  Standard  Specific  Angle  of 

Deviation  (dq)  Gravity  (s)  Repose  (<t>) 


Glass  Beads  0.220 
Ottawa  Sand  0.585 


1.08 

1.16 


2.49 

2.62 


22° 

32.5° 


Only  for  a/rQ  =  1  does  an  analytical  solution  exist  for  relating  rate  of 
flow  tc  the  coefficient  of  permeability  and  .  The  variation  of 

the  vertical  velocity  component,  V  ,  along  the  bed  is,  from  equation  (6) 


V  » « “ 

Y“0  ix 


(?) 


in  which  Ah  =  h_  -  1c  .  The  total  flow  rate,  Q  is  obtained  by  inte¬ 
grating  equation  (7)  over  the  area  of  the  interface 


Q  = 


y=0 


dA  =  k  fcAh 

it 


dx 

x 


in  which  b  is  the  width  of  flow  passage.  Upon  integrating  and  reducing 


Q  =  0.512  bK  Ah 


(8) 


The  test  procedure  consisted  of  setting  the  flow  rate  Q,  waiting  until 
a  state  of  equilibrium  was  reached,  and  then  measuring  all  pertinent 
quantities.  The  discharge  Q  was  measured  gravimetrically,  the  change  in 
piezometric  head  Ah  was  read  off  the  scale  in  Figure  5  and  the  coefficient 
of  permeability  was  computed  from  equation  (8).  At  60°F  the  coefficient 
of  permeability,  k,  is  0.00103  ft/sec  and  0.00465  ft/sec  for  the  glass 
beads  and  Ottawa  sand,  respectively. 


Boiling  and  Heaving 

Tests  were  conducted  in  order  to  study  boiling  and  heaving  near 
point  B  on  the  cylinder.  Of  particular  interest  is  the  behavior  of  the 
flow  near  singular  point  B  for  the  case  in  which  the  cylinder  is 
partial  ly-  imbedded . 

Cylinder  Half-Imbedded 


With  the  cylinder  half-imbedded  in  the  bed  material  the  sequence 
of  events  from  initial  boiling  to  piping  were  observed.  The  glass  beads 
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and  the  Ottawa  sand  were  used  as  bed  materials.  The  bed  was  screeded 
and  made  plane  as  shown  in  Figure  5.  The  test  procedure  consisted  of 
increasing  the  piezometric  head  In  2  slowly  until  boiling  and,  finally, 
piping  occurred.  As  the  water  level  h0  was  raised  the  first  sand  move¬ 
ment  at  B  was  localized  boiling  for  the  glass  beads  and  slow  heaving  for 
the  Ottawa  sand.  For  the  glass  beads,  the  boiling  was  localized  at 
first  but  then  nearly  uniformly  distributed  over  the  bed  adjacent  to 
the  cylinder  at  B.  As  the  piezometric  head  was  slowly  increased  the 
boiling  became  violent  to  the  extent  that  the  glass  beads  began  to  move 
in  mass,  forming  a  mound  adjacent  to  the  cylinder.  For  the  Ottawa  sand 
the  mound  grew  from  initial  heaving  rather  than  boiling  but  the  shape  of 
the  mound  is  similar  to  that  for  the  glass  beads.  A  photograph  of  the 
mound  at  B  and  the  depression  at  D  for  the  glass  beads  is  shown  in 
Figure  6a.  The  mound  for  the  Ottawa  sand  is  shown  in  Figure  6b.  Each 
mound  grows  in  size  but  is  similar  in  shape  until  a  pipe  forms. 

As  stated  above,  piping  does  not  immediately  follow  boiling. 

There  is  definitely  a  self-correcting  mechanism  present,  up  to  a  limit, 
during  the  heaving  process.  For  several  mound  sizes  no  piping  will  occur 
provided  Ah  is  maintained  constant.  The  limiting  case  occurs  when  a 
pipe  develops.  Hence  there  is  a  factor  of  safety  (S.F.)  in  designing  for 
a  critical  gradient.  This  factor  varies  from  1.16  to  1.36,  as  shown  in 
the  table  below. 

The  angle  that  the  sloping  portion  of  the  mound  makes  with  the 
horizontal  can  be  explained  qualitatively.  This  angle  is  always  less 
than  the  natural  angle  of  repose  of  the  particular  bed  material  because 
of  the  water  flowing  out  of  the  bed.  The  value  of  this  angle  before  piping 
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occurs  is  15°  for  the  glass  beads  and  20°  for  the  Ottawa  sand.  As 
discussed  by  Haefeli  (3),  a  bed  is  rendered  less  stable  with  flow  out 
of  the  bed  and  more  stable  with  the  flow  into  the  bed.  Such  an  effect 


causes  the  bed  slope  to  be  less  than  the  natural  angle  of  repose  when 
fiow  is  out  of  the  bed  and  more  than  the  natural  angle  of  repose  when 
flow  is  into  the  bed.  Both  of  these  cases  are  borne  out  by  the  mild 
slope  on  the  mound  and  the  steep  slope  on  the  depression  near  D.  The 
angle  that  the  depression  makes  with  the  horizontal  is  nearly  34°  for  the 
glass  beads  and  42°  for  the  Ottawa  sand. 

Cylinder  Partially-Imbedded 

The  main  objective  for  the  tests  with  the  cylinder  partially- 
imbedded  was  to  observe  the  behavior  of  the  bed  in  the  vicinity  of 
singular  point  B.  Tests  were  run  with  a/rQ  =  0.25,  0.50,  and  0.75. 

Tests  with  a/rQ  =  0.25  and  0.50  were  unsatisfactory  as  piping  occurred 
before  a  mound  formed  completely  across  the  channel.  All  data  presented 
are  for  a/rQ  =  0.75.  Both  the  glass  beads  and  the  Ottawa  sand  were  used 
as  bed  materials.  The  bed  was  screeded  plane  at  a  level  l/2  in  below 
the  centerline  level  of  the  cylinder,  Figure  7.  The  bed  adjacent  to 
singular  point  B  was  observed  as  the  piezometric  head  h^  was  increased 
slowly. 

Boiling  was  net  observed  immediately  as  one  might  expect  from  the 
theoretical  solution.  In  fact,  for  the  Ottawa  sand,  initial  movement 
was  by  heaving,  as  with  a/rQ  =  1.  Boiling  or  heaving  do  occur,  of  course, 
at  a  value  of  Ah  less  than  that  required  for  the  half-imbedded  cylinder. 
As  the  water  level  h^  increased  the  boiling  for  the  glass  beads  became 
more  violent  until  a  mound  was  formed.  The  interesting  feature  is  that 
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the  mound  is  curved,  unlike  the  rectilinear  features  of  the  mound  for  the 
half-imbedded  cylinder.  In  Figure  8  are  shown  photographs  of  a  typical 
mound  for  the  glass  beads  and  for  the  Ottawa  sand.  Upon  close  scrutiny 
it  is  noticed  that  the  mounds  intersect  the  cylinder  at  essentially 
right  angles.  Nature  apparently  corrects  for  the  anomaly  of  non¬ 
orthogonality  at  such  singular  points  as  B.  With  a  further  increase  in 
Ah,  piping  finally  occurs.  Values  of  Ah  for  the  initial  boiling  and 
piping  and  values  of  S.F.  are  given  in  the  table  below. 

The  self-correcting  mechanism  is  also  present  for  a/r^  <  1.  The 
initial  formation  of  a  mound  takes  place  fairly  rapidly  once  boiling 
commences  but,  if  h 2  is  held  constant  the  heaving  of  the  mound  ceases. 
Only  if  h2  is  increased  further  does  any  movement  occur.  Hence  the 
mound  is  stable,  up  to  a  limit,  for  a  constant  h2-  The  heaving  process 

for  a/r  <  1  is  similar  to  that  for  a/r  =  1. 

o  o 

Piping  and  Scour 

As  the  piezometric  head,  h2>  is  further  increased  the  mound  and 
depression  grow  in  size  until  a  pipe  is  formed.  Once  a  pipe  develops 
under  the  cylinder  an  unstable  condition  exists  such  that  scour  is  immi¬ 
nent.  The  piping  and  scour  for  the  partially-imbedded  cylinder  is  so 
similar  to  that  for  the  half-imbedded  one  that  no  distinction  between  the 
two  is  necessary. 

Motion  pictures  were  taken  of  the  entire  piping  and  scour  sequence. 
Once  piping  appeared  imminent  to  the  observers  a  stop-clock  was  started 
and  the  inlet  flow  was  stopped.  The  motion-picture  camera  was  started 
and  the  entire  sequence  photographed.  The  sequence  will  be  described 


a.  Glass  Beads  as  Bed  Material 


Figure  8.  Profile  of  Mound  for  a/r 
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with  reference  to  Figures  9  a,  b,  c,  and  d.  In  Figures  6  and  8  the 
mound  and  depression  are  in  a  stable  state.  In  Figure  9a  a  pipe  is 
seen  adjacent  to  the  cylinder  at  the  near  wail.  At  this  instance  scour 
is  inevitable  as  a  channel  is  forming  under  the  cylinder.  In  Figures  9b 
and  c  is  shown  the  development  of  the  scour  channel  beneath  the  cylinder. 
Since  the  flow  into  the  apparatus  is  stopped  just  before  piping,  and 
the  piezometric  head  h ^  falls  rapidly,  the  velocity  out  of  the  channel 
decreases  to  a  value  below  which  scour  will  occur.  At  this  condition 
the  bed  material  can  no  longer  be  carried  out  of  the  hole  and  sloughing 
back  into  the  channel  occurs.  The  sloughing  continues  until,  finally, 
with  h^  =  h^,  the  slopes  of  the  bed  on  the  left-  and  right-hand  sides 
reach  the  natural  angle  of  repose  of  the  bed  material.  This  final  state 
is  shown  in  Figure  9d. 


Conclusions 

Singular  points  which  occur  at  imposed  conditions  of  nonortho¬ 
gonality  of  flow  lines  and  potential  lines  are  not  as  critical  as 
mathematical  analysis  indicates.  Boiling  does  not  occur  immediately 
but,  when  it  does,  the  bed  takes  on  a  shape  such  that  orthogonality  is 
satisfied  at  these  singular  points.  A  safety  factor  defined  by  the 
ratio  of  the  piezometric-head  gradient  required  to  cause  piping  to  that 
required  to  cause  boiling  varies  from  1.16  to  1.57  for  all  tests 
involving  a  cylinder  being  either  partially-  or  half-imbedded. 
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Table.  Results  for  Initial  Boiling  and  Piping 


Ah  (in) 


a/ro 

Bed 

Material 

Initial 

Boiling 

Initial 

Piping 

S.F. 

1.00 

Glass  Beads 

6.34 

7.34 

1.16 

1.00 

Glass  Beads 

5.72 

7.59 

1.32 

1.00 

Glass  Beads 

6.50 

7.70 

1.18 

1.00 

Glass  Beads 

5.90 

7.49 

1.27 

1.00 

Glass  Beads 

5.82 

7.35 

1.26 

1.00 

Ottawa  Sand 

8.52 

10.27 

1.21 

1.00 

Ottawa  Sand 

7.53 

10.25 

1.36 

0.75 

Glass  Beads 

3.75 

5.78 

1.54 

0.75 

Glass  Beads 

3.84 

6.02 

1.57 

0.75 

Ottawa  Sana 

6.53 

8.03 

1.23 
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